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In 1974, Gehring posed the problem of minimizing the length of two linked curves 
separated by unit distance. This constraint can be viewed as a measure of thickness 
for links, and the ratio of length over thickness as the ropelength. In this paper 
we refine Gehring's problem to deal with links in a fixed link-homotopy class: we 
prove ropelength minimizers exist and introduce a theory of ropelength criticality. 

Our balance criterion is a set of necessary and sufficient conditions for criticality, 
based on a strengthened, infinite-dimensional version (Theorem 5.4) of the Kuhn- 
Tucker theorem. We use this to prove that every critical link is C' with finite 
total curvature. The balance criterion also allows us to explicitly describe critical 
configurations (and presumed minimizers) for many links including the Borromean 
rings. We also exhibit a surprising critical configuration for two clasped ropes: 
near their tips the curvature is unbounded and a small gap appears between the two 
components. These examples reveal the depth and richness hidden in Gehring's 
problem and our natural extension. 

57M25; 49Q10, 53A04 



1 Introduction 

Suppose that A and B are disjoint linked Jordan curves in 

which lie at a distance 1 from each other. 
Show that the length of A is at least In. 

—Fred Gehring, 1974 

Gehring's problem, which appeared in a conference proceedings [7], was soon solved 
by Marvin Ortel. Because Oriel's elegant solution was never published, we reproduce 
it here with his permission: Fix any point a G A; the cone on A from a is a disk 
spanning A . Since A and B are linked, B meets this disk at some point b £ B, lying 
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on a chord of A. Because Dist(A, b)> \ , projecting A to the unit sphere 5 around b 
does not increase its length. The projection is a closed curve joining two antipodal 
points on S, and so has length at least It: . (Further proofs and generalizations to linked 
spheres in higher dimensions were published by Edelstein and Schwarz [11], Osserman 
[28] and Gage [12, 13].) 

The unique minimizing configuration for Gehring's problem is a Hopf link consisting of 
two congruent circles in perpendicular planes, each passing through the other's center. 
This leads to a natural question: what are the length-minimizing shapes of other link 
types when the different components stay unit distance apart? This constraint prevents 
different components from crossing each other, but we cannot expect to fix the link 
type exactly. Instead, the natural setting for this problem is Milnor's notion of link 
homotopy: two links are link-homo topic if one can be deformed into the other while 
keeping different components disjoint. Clearly one link can be deformed into another 
while keeping all components at unit distance if and only if they are link homotopic. 

We will define the link-thickness of a link to be the minimum distance between different 
components. The problem we consider is then to minimize length in a link-homotopy 
class, subject to the constraint of fixed link-thickness. Equivalently, we could minimize 
the link-ropelength of the link, meaning the quotient of length over thickness. 

In [4], we found length-minimizing links under a similar constraint: that a normal tube 
of diameter one around the link stay embedded. It is easy to see that the examples 
constructed there (like the one in Figure 4) are also global minima (in their respective link- 
homotopy classes) for the Gehring problem. The focus of this paper will be on critical 
configurations. Our main result is a balance criterion (Theorem 6.1, Corollary 6.3), 
which states that a link is link-ropelength critical if and only if the tension force in the 
curve is balanced by a system of compressive forces between pairs of points on different 
components of L realizing the minimum distance. 

This balance criterion is based mainly on an improved, infinite-dimensional version 
(Theorem 5.4) of the Kuhn-Tucker theorem on constrained optimization, which is 
essentially a very general method of Lagrange multipliers. The other key technical 
element is a careful application of Clarke's differentiation theorem for min-functions 
(Theorem 3.1). 

The direct method shows that there is a (rectifiable) minimizer for link-ropelength in 
each link-homotopy class. An interesting problem is to determine the regularity of these 
minimizers or other critical points. The previously known minimizers were C^'^ but 
not . Our balance criterion allows us to prove that all link-ropelength-critical curves 
are with finite total curvature (Proposition 6.5). 
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We next consider generalized links, which may include open components with con- 
strained endpoints, or which avoid fixed obstacles. After extending our balance criterion 
and existence results to this setting, we analyze the problem of the simple clasp. A 
clasp consists of two linked arcs whose endpoints are constrained to parallel planes 
(as in Figure 10). A generalization to clasps of different opening angles provides a 
model for the strands of rope in a woven cloth or net. The balance criterion lets us 
construct explicit critical configurations (Theorem 9.5) of these generalized links; we 
conjecture they are the length-minimizers subject to the constraint that the arcs remain 
at unit distance from each other. Our critical clasp has a number of surprising features, 
including a point of infinite curvature and a small gap (at the center of the clasp) between 
the tubes around the two components. This configuration is C^'^/^ and may represent 
the worst regularity of any critical curve. 

We end by constructing a ropelength-critical configuration (and presumed minimizer) 
for the Borromean rings. In all the other known critical configurations for closed links, 
each component is a convex plane curve built from straight segments and arcs of circles. 
In our Borromean rings, the components are still planar, but are nonconvex and are 
built from different pieces including parts of a clasp curve. In a sense, this is the first 
nontrivial example of a ropelength-critical link. 

Our methods will have a number of other applications. In particular, we have used 
them to describe critical configurations for the "standard" ropelength problem for knots 
and links: minimize the length of a link subject to the constraint that the normal 
neighborhood of unit diameter remains embedded. We will publish these results in a 
sequel [2] to the current paper. We can also consider minimization not of length but 
of other objective functions like elastic bending energy, again subject to a thickness 
constraint. Analogs of our balance criterion may be useful in describing other flexible 
mechanisms, such as thick surfaces. 

We note that von der Mosel and Schuricht [33] have used a similar approach (via 
Clarke's theorem and a functional-analytic version of Lagrange multipliers) to derive 
necessary, but not sufficient, conditions for criticality for the ropelength functional 
of [4]. We will treat the same functional in our forthcoming sequel [2] and will offer a 
comparison of the two methods there. We also note that Starostin has given [34] an 
independent derivation of the tight clasp of Section 9, though he does not prove that it 
is critical. 
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2 Link-thickness for closed links 

In order to reformulate Gehring's problem, we first establish some basic terminology. 
Remember that a compact, oriented 1 -manifold- with-boundary M is a finite union of 
components, each of which is homeomorphic to a circle §^ or an interval [0, 1] . 

Definition A parametrized curve is a mapping from a compact, oriented 1 -manifold- 
with-boundary M to . Two parametrized curves are equivalent if they differ by an 
orientation-preserving reparametrization (thatis, by composition with an orientation- 
preserving self-homeomorphism of M). A curve L in is an equivalence class of 
parametrized curves. We say L is closed when each component of its domain M is a 
circle, that is, when its boundary dL is empty. 

Even though our curves may have self-intersections, we will usually refer to points on 
the curve as if they were simply points of its image in R^ . The meaning should be clear 
from context. 

The length Len(L) of any curve L is defined to be the supremal length of all polygons 
inscribed in L. A curve has finite length, or is rectifiable, if and only if it has a Lipschitz 
(that is, C"'^ ) parametrization. One such parametrization is then by arclength s. Any 
rectifiable curve has a well-defined unit tangent vector T = AL/As almost everywhere. 

Definition The link-thickness LThi(L) of a curve L is the minimum distance between 
points on different components of L. This is the supremal e for which the {^/i}- 
neighborhoods of the components of L are disjoint. 

For now, we will consider only the case of closed curves, where each component is 
a circle. (We will deal with generalized links — with endpoint constraints — later in 
Section 8.) So suppose we start with a closed curve L and we want to minimize length 
under the constraint that the link-thickness remains at least one. Since we can rescale 
any link to have LThi > 1 , this problem is the same as minimizing (link-) ropelength, 
the quotient of length by link-thickness. 

The thickness constraint naturally prevents different components from passing through 
each other, but does not prevent any given component from changing its knot type 
through self-intersections. This is the setting for Milnor's work on link homotopy: 

Definition A link is a closed curve with disjoint components. The link-homotopy class 
of a link L, denoted [[L]] , is the set of curves homotopic to L through configurations 
that keep different components of L disjoint. 
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Note that, for our purposes, configurations of L where some components have self- 
intersections are still considered to be links, and are included in [[L]] . 

For two-component links, Milnor [25] showed that linking number is the only link- 
homotopy invariant. For links of many components, the topological situation is more 
complicated, but a complete classification of links up to link homotopy was provided 
by Habegger and Lin [16]. We will prove in Section 6 that in every link-homotopy 
class there is a curve minimizing ropelength. We show these minimizers are always , 
though our examples suggest that they may not always have bounded curvature. 

3 The derivative of link-thickness 

We want to define critical configurations of L subject to the thickness constraint 
LThi(L) > 1 . Because LThi is defined as the minimum of a collection of distances 
between points on different components, the equation LThi > 1 acts like a collection 
of many constraints. To make this notion precise, we will apply a theorem of Clarke to 
compute the derivative of LThi as we vary the curve L. 

Given any curve L, let L*^^ be the compact set of all unordered pairs {x, y} of points on 
distinct components of L. The link-thickness of L is simply the minimum over L^^^ of 
the distance function Dist{x,y} := |y — x| . 

We often want to consider a continuous deformation of a curve L: fixing any 
parametrization / of L, that means a continuous family of parametrized curves with 
/o =/. (When we reparametrize L, we apply the same reparametrization to Lt at all 
times t.) We assume that Lt is in t; the initial velocity of L, will then be given by 
some (continuous, -valued) vector field ^ along L. We let VF(L) denote the space 
of all such vector fields. Formally, these are sections of the bundle /*rR^ pulled back 
from the tangent bundle of ^ by the parametrization / of L. Identifying any tangent 
space to with M'' itself, this is simply a map from the domain M to . Again, 
when we reparametrize a curve L, we apply the same reparametrization to any vector 
field ^. 

Consider a curve L with LThi(L) > 0. If is a continuous deformation of L, with 
initial velocity given by some ^ G VF(L), then for each pair {x, v} G L^^\ we clearly 
have 



(Even if L is not embedded, the condition LThi(L) > implies x and y cannot coincide 
in M.^ , so this formula is always meaningful.) 



6^Dist{x,y} := ;r;\y - x\ 



\y-x\ 
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A function like LThi, defined as tlie minimum of a compact family of smooth functions, 
is sometimes called a min-function. Clarke's differentiation theorem for min-functions 
says that — ^just as in the case when the compact family is finite — the derivative of 
a min-function is the smallest derivative of those smooth functions that achieve the 
minimum. More precisely, specializing [6, Theorem 2.1] to the case we need, we have: 

Theorem 3.1 (Clarke) Suppose for some compact space K and some e > 0, we 
have a family of functions fk : (— e, e) — > M, for k ^ K. Suppose further thatfk{t) 
and flit) are lower semicontinuous on K x (—e,e). Letf(t) := mirik^zK fk(t)- Thenf 
has one-sided derivatives, and 

df 



= min//(0), 

t=0 keKo 



where Kq := {k £ K : fk(0) = /(O)} is the set of k for which the minimum in the 
definition off is achieved when t = 0. 

To apply this theorem to thickness, suppose we have a variation Lt of the curve L, 
and let ^ G VF(L) be its initial velocity. The link-thickness LThi(L) is written as a 
minimum over K = L(2) of the pairwise distance. Clarke's theorem picks out those 
pairs achieving the minimum: Kq is the set of pairs achieving the minimum distance 
LThi(L). 

In rigidity theory, the vertices of a tensegrity framework are joined by bars whose length 
is fixed, cables whose length can shrink but not grow, and struts whose length can grow 
but not shrink (compare Roth and Whiteley [31]). Thus, we borrow the term "strut" to 
describe the pairs in A'o : 

Definition An unordered pair of points {x, y} on different components of L is a strut 
if Ij — Jc| = LThi(L). The space of all struts of L is denoted Sti-ut(L) C L^^K 

Struts correspond to points of contact between tubes around the different components 
of L. Our balance criterion will show how the segment xy can be viewed as carrying a 
force pushing outward on its endpoints. 

Applying Clarke's theorem to link-thickness, we get: 

Corollary 3.2 For any curve L, and any variation vector field ^ G VF(L), the 
(one-sided) first variation of link-thickness is 

5t LThi(L) = min 6cDist{x, y| . 

? Stmt(L) ^ 



Qeometry & Topology 10 (2006) 



Criticality for the Gehring link problem 



2061 



Note that (5^Dist{x, y} is a continuous function of x and y, and for any fixed {x, y} is a 
linear function of the variation ^, being the derivative of a smooth function. 

Therefore we can collect these into a linear operator As = (5Dist from VF(L) to the 
space C(Strut(L)) of continuous functions on struts, defined by 

(AsOCl^,)'}) := <5pist{x,3;} = {^y -^x,y-x). 

Borrowing again from rigidity theory, where the analogous A is called the rigidity 
matrix, we will call As the rigidity operator for link-thickness. 

The corollary above can be rephrased to conclude that a variation ^ decreases LThi(L) 
to first order if and only if As^ takes at least one negative value on Strut(L). 

Note that, while the corollary says that link-thickness has a directional derivative 
6^ LThi in each direction ^, the operator 5'^ LThi is not linear in ^. For instance, when 
one component of a Mnk is between two others, it is easy to have both S'^ LThi < and 
S^^ LThi < 0. We write the superscript to emphasize that these are only one-sided 
derivatives. There is, however, a form of superlinearity: 

Corollary 3.3 For any curve L and any ^, G VF(L), we have 
(5+^^ LThi(L) > S+ LThi(L) + S+ LThi(L). 

Proof This follows immediately from the linearity of As and the general fact that 
+ g) > niin/ -I- ming. We have 

minAs(C + v){x, y} = min(AsC + Asri){x, y} 

> rmnAs^{x, y} + min As?7{x, y}, 
where the minima are taken over all {x, y} G Strut(L) . □ 

We will be interested in the adjoint Ag of the rigidity operator, so we first consider the 
dual function spaces. By the Riesz representation theorem, we know that C* (Strut) is 
the space of signed Radon measures on the space Strut(L) of struts. Similarly VF*(L) 
is the space of what we will call forces along L, namely vector-valued Radon measures 
on L. 

The adjoint operator A^ now associates to any measure fi on struts a force A^/i along L. 
Geometrically, each pair {x, y} acts along the chord xy, outward at each of its endpoints. 
In formulas, 

/ CdAli^ = [ AsCd/x 

Jl i Strut 

= / / Uy,r^)^f^^^y^^ 
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where we have hfted /^i to a symmetric measure ^^{x,y) on ordered pairs. Physically, 
we think of ji as giving the strengths of compressive forces within the struts, and as 
the operation that integrates these strut forces to give a net force along the curve L. 

4 First variation of length and finite total curvature 

The objective functional we consider in this paper is simply the length Len(L) of a 
curve. Since our curves might not be smooth, we need to carefully examine the first 
variation of length. 

Let L be a rectifiable curve parametrized by arclength s, with unit tangent vector T . 
Suppose Lt is a variation of L under which the motion of each point x G L is smooth in 
time with initial velocity , and ^ € VF(L) is a Lipschitz function of arclength. Then 
the standard first- variation calculation shows that 



where ^' = d^/ds is the arclength derivative, defined almost everywhere along L. 
If L is smooth enough, we can integrate this by parts to get 



(In the boundary term, the sign is chosen to make ±T point inward at x.) In fact, not 
much smoothness is required: as long as T is a function of bounded variation, we can 
interpret T' as a measure, and the formula holds in a sense we will now explore. 

Following Milnor [24], we recall that the total curvature of a polygon is just the sum of 
its (exterior) turning angles, and we define the total curvature of any curve to be the 
supremal total curvature over all inscribed polygons. A rectifiable curve L has finite 
total curvature if and only if the unit tangent vector T = L'{s) is a function of bounded 
variation. Sometimes the space of all such curves is called W^'^^ or BV^ but we will 
call it FTC. (See Sullivan [36] for a survey of results on FTC curves.) 

If L G FTC , it follows that at every point of L there are well-defined left and right 
tangent vectors T± ; these are equal and opposite except at countably many points, the 
corners of L. (See, for instance, Royden [32, Sect. 5.2].) 

If L is FTC , its tangent T has a distributional derivative /C with respect to arclength: a 
force (an -valued Radon measure) along L that we call the curvature force. 
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The curvature force has an atom (a point mass or Dirac deUa) at each corner x e L, with 
lC{x} = + r_(x). On a arc of L, the curvature force is /C = dJ = kN As 

and this is absolutely continuous with respect to the arclength or Hausdorff measure 
d5 = 

When L has boundary, we choose to include in K, an atom at each endpoint of L, with 
mass 1 and pointing in the inward tangent direction. This means we need no boundary 
terms in the formula (5^Len(L) = — (C) d/C) . 

We say that a vector field ^ along L is smooth if is a smooth function of arclength. 
(The arclength parametrization of any rectifiable curve is essentially unique, so this 
makes sense.) The set of all smooth vector fields will be denoted VF°°(L). 

The first variation JLen(L) can be viewed as a linear functional on smooth vector 
fields ^ G VF°^(L). As such a distribution, it has order zero, by definition, if 
5gLen(L) = /^(T, d^ is bounded by C sup^ |^| for some constant C. This happens 
exactly when we can perform the integration by parts. 

We collect these results as: 

Lemma 4.1 Given any rectifiable curve L, ttie following conditions are equivalent: 

(a) Lis FTC. 

(b) The first variation JLen(L) has distributional order zero. 

(c) There exists a curvature force measure K, along L such that d^Len{L) = 

-/^(e,d/c). 

An FTC curve L is exactly when it has no comers, that is, when JC has no atoms 
(except at the endpoints). It is furthermore C^'^ when T is Lipschitz, or equivalently 
when JC is absolutely continuous (with respect to arclength) and has bounded Radon- 
Nikodym derivative dlC/dTL^ = kN. In previous work on ropelength (for example 
[4, 14]), the thickness measure had an upper bound for curvature built in, meaning 
that any curve of positive thickness was automatically C''^ . This is not true for the 
link-thickness, so we do not expect the same regularity results to hold here. 

5 Constrained criticality and the Kuhn-Tucker theorem 

We will review constrained minimization problems in a finite setting, before generalizing 
to the setting we will need for our ropelength problems. Suppose we want to minimize 
a function / : M" ^ M inside the admissible region defined by a finite collection 
of inequality constraints > 0. A constraint gi is active at p £ W if gi{p) = 0. 
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Definition We say that /? is a constrained critical point for minimizing / if, for any 
tangent vector v at p with D^f < 0, we have Dygi < for some active gi. That is, p is 
critical if there is no direction v € M" that reduces / to first order while preserving all 
constraints to first order. 

Note that the criticality conditions for minimizing / and — / are quite different; in 
particular local maxima for / are rarely critical points for minimizing/, while local 
minima for / usually — though not always — are. 

Example 5.1 Suppose we minimize f{x, z) := x on the halfplane x > in , subject 
to 

gi :=(x2-l)3-z>0, g2:=z>0. 

The admissible region has an outward-pointing cusp, shown in Figure 1. The tip of this 
cusp, at p = ( 1 , 0) , is the global minimum of / over the admissible region, but it is 
not critical: the directional derivatives in the direction v = (— 1, 0) are D,/ = — 1 but 
Dvgi = 0. 




Figure 1: In this illustration of Example 5.1, the admissible region for the two constraints 
gi = (x^ — 1)^ — z > and g2 — z > is shaded. The Mangasarian-Fromovitz constraint 
qualification fails at the cusp point p ~ (1,0) because Vgi and V^2 are equal and opposite 
there. 

To deduce that a local minimum of / is critical according to our definition, an 
additional regularity hypothesis will be required. However, critical points can be exactly 
characterized by a Lagrange multiplier theorem (compare [17]): 

Theorem 5.2 (Modified Kuhn-Tucker Theorem) A point p is constrained-critical 
for minimizing f if and only if the gradient V/ is a positive linear combination of the 
gradients Vgj of the constraints gj active at p. 
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The geometric intuition behind this theorem is easy to understand: Only the active 
constraints matter, and being inequahty constraints they can only act positively. If there 
were some component of —V/ not canceled by the Vgj, that would give an admissible 
direction to move which decreases/. 

Unlike in the classical Kuhn-Tucker theorem, we do not need additional regularity 
hypotheses on the point p , which may surprise those familiar with optimization theory. 
The explanation is that we are interested in critical points, while the classical theorem 
deals with minima of/. And as we saw above, not every minimum of / is critical. But 
just as in the classical theory, criticality will be guaranteed if we add the hypothesis that 
the Mangasarian-Fromovitz constraint qualification [20] holds for a local minimum. 

Definition A point p is constraint-qualified (in the sense of Mangasarian and Fro- 
movitz) if there is a direction v such that for all constraints gj active at p we have 
D,gj > 0. 

We note that this condition fails at the point /? = (1, 0) in Example 5.1 above, which 
was minimal but not critical. 

Proposition 5.3 If p is a local minimum for/ when constrained by {gi > 0} , and p 
is constraint-qualified, then p is constrained-critical for minimizing f . 

We have omitted proofs of the theorem and proposition above because they are standard 
and are also special cases of our infinite-dimensional generalizations below. 

5.1 A generalized Kuhn-Tucker theorem 

Note that in Theorem 5.2, the functions / and gi might as well be replaced by linear 
functions — their differentials at p. We view this as the linear- algebraic core of the 
Kuhn-Tucker theorem. 

We will now derive an infinite-dimensional version, where the linear functional / is 
defined on an arbitrary vector space X, and the finite family of constraints gi is replaced 
by a family Ay, where y ranges over some compact space Y. 

While our theorem does not mention optimization directly, it will be the engine that 
drives all of the optimization theorems of this paper. 

As usual, we let C{Y) be the Banach space of continuous functions on Y with the sup 
norm || • || , and let P C C{Y) be the closed positive orthant consisting of nonnegative 
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functions. Then the dual space C*(F) consists of all signed Radon measures on Y , and 
P* c C*(F) is the cone of positive measures. 

Note that any function z S C(F) can be decomposed into positive and negative parts: 
z = z+ — z~ with € P- Then we have \\zr \\ = Dist(z, P). 

Theorem 5.4 Let X he any vector space and Y he a compact topological space. For 
any linear functional f on X and any linear map A : X — > C{Y), the following are 
equivalent: 

(a) There exists e > such that ||(A0"|| > £ for all ^ € X withfiO = -1. 

(h) There exists a positive Radon measure n e P* such that /(O = /i(AO for all 

This theorem is comparable to the generalized Kuhn-Tucker theorem of Luenberger 
[19, page 249]. His theorem, restated to apply to the linear Gateaux differentials (f 
and A) of the original objective and constraint functions on X, says: 

Theorem 5.5 (Luenberger) Let X and Z he vector spaces, with a norm given on Z, 
and let P C Z he a closed convex cone with nonempty interior. Letf : X ^ M. he a 
linear functional and A : X — > Z be a linear map. Assume that whenever A^ £ P we 
havefiO ^ 0' f^^f ^^^s in the interior of P for some G X. Then there exists 
HeP* such thatfiO = KM) for all ^ G X. 

While our version applies only to the case Z = C{Y), our hypotheses (a) on / and A 
are somewhat weaker than those imposed by Luenberger — they are necessary as well as 
sufficient for (b) the existence of /i . 

To understand our overall strategy, consider the linear map (f, A) : X ^ M x C(F). As 
we will see below, (a) implies that the image of (/, A) avoids the interior of the orthant 

R- X P. 

To gain some intuition, let us specialize to the case where X = W" . We can rephrase (b) 
to say that some vector in the kernel of the adjoint map (f, A)* is in xP. When Y has 
finite cardinality n, we put the standard Euclidean inner product on M x C(F) = M"^^ , 
and identify this space with its dual. Then the kernel of (f,A)* and the image of (f,A) 
are orthogonal complements in M"+^ . The standard Farkas alternative (see Figure 2) 
says that, given any closed orthant in M"+^ , it must intersect one out of any pair of 
orthogonal complements. Our argument in the general case, when F might be infinite 
and X infinite dimensional, will be guided by this intuition. 
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Figure 2: One version of the Farkas alternative states that, given any closed orthant O in a inner 
product space, it must intersect at least one out of any pair of orthogonal complements / and . 

Proof of Theorem 5.4 One direction is easy: suppose we have a positive Radon 
measure /x so that for each ^ G X, 

/(o = ^Aed/i. 

For any ^ G X with/(0 = -1, write z := G C{Y). We have / zd// = -1, and 
since ^ is a positive measure, we can replace the function z with its negative part to 
conclude that /z^ d^u > 1. Furthermore /U has finite mass mass(//) := / d/n < oo by 
the Riesz theorem. Therefore 

Dist(z,P) = llz^ll > l/mass(//) > 0. 

This completes the proof that (b) implies (a). 

To prove the converse, first give M x C(F) the Euclidean combination of the sup norms 
on M and C(y): 

\\{.a,g)\\ = ^/a^+\W. 

Now consider the orthant O := [—1, oo) x P. Our hypothesis (a) implies that there 
is positive distance between O and the image / := {f,A){X) of the linear map 
(f,A). Take sequences (/(^Oi^CO in I ^nd {ti,Zi) in O, whose pairwise distance 
approaches Dist(/, O) . That is, setting 

v;:= {ti-mi),Zi-Ai,) 

we have ||v, || ^Dist(/, O). 

We first claim that we can assume that v,- G x P. Certainly we can assume the 
equality Zi = , since this positive part of the function A^,- realizes the distance 

||z,- - A^-ll = ||(A6-)"|| = Dist(Ae-,P). Then z; - A^- = -(A^)" G P- Similarly, 
Dist(/'(,^,), [— 1, oo)) is —1 —fiii) if this is positive (and is zero otherwise), so we may 
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assume ti = min(— Thus f,- <f{ii), so ti —f(£,i) < 0. This proves the first 
claim. 

We now have a geometric problem: from Figure 2 we see there is a special case where 
both the image of (/", A) and its orthogonal complement lie on the boundary of x P. 
If this happens, then the closure I intersects either the subspace M x {0} or {0} x P. 
The second case does not trouble us, but the first would cause us problems later; we 
now show that our assumption (a) rules it out. To do so, we think about the setup above 
geometrically: if I intersects M x {0}, then we expect that ti — ^ 0- 

Thus our second claim is that we can assume the U —fiCi) are uniformly negative. 
If not, lim/(^,) < — 1 , so without loss of generality, we can rescale ^, down so that 
fi^i) = — 1 . That means (— 1, A^,) S /. By hypothesis (a) we know 

di := Dist((-l,Ae), O) = Dist(Ae-, P) > e 

for some fixed e > 0. Since we are using the Euclidean combination of the norms 
on M and C(F), the distance from any rescaling by k of (— 1, A^,) to O is given by the 
Euclidean distance from {—k, k\\A^f ||) to (—1,0). And we can use plane geometry to 
see that rescaling by 1 — brings us closer to O : 

Dist((l - e2)(-l,A^,), O) = di^l -le^ + e'^{\ + l/d}) < di^/l -e^ + e"^. 

We can always assume that e < 1, so the constant Vl — + e'^ is less than 1. 
Therefore 

Dist(/,0) < limDist((l -e2)(-l,A6), O) 

< limDist((-l,A^,), O) 
= Dist(/,0). 

This contradiction proves the second claim. 

We have proved that the v, are in x P. Using the Hahn-Banach theorem, for 
each / we can find a linear functional (c;, Ui) G M x C*(F) that vanishes on /, satisfies 
{ci, fi)(v,) = 1, and has norm 

||(c;,z.O|| = 1/Dist(v,- + /, (0,0)). 

Because the M-components of v; are uniformly negative, so are the c,-. 

Using Alaoglu's theorem, the (c,-, ui) have a subsequence converging in the weak* 
topology to a limit functional (c, v) ; we have c < and its norm is bounded above by 
1 / lim Dist(vi, -/) = 1 /Dist(/, O) . 
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Setting ^ := i'/\c\ G C*(F), we claim this will be the Radon measure in statement (b). 
By construction, (— 1, vanishes on /, meaning that for G X, we have 

(Notice that we have used the additional geometric information that / does not approach 
M X {0} in an essential way; if it did, then c would vanish, and we could not rescale v 
by l/|c| to obtain the equation above.) 

It remains only to show that ji is positive. In an inner product space, this would be 
obvious: each Ui would be positive (since it was dual to a positive Zi — A^, ), and u 
would be a limit of positive measures. But our Vj were constructed implicitly by the 
Hahn-Banach theorem, and so might include negative pieces. We now address this 
problem. 

We can decompose each Vi into its positive and negative parts Vj = vf' — , with 
mass(t',) = mass(f,^) + mass(i/,~). In order to show v is positive, we will prove that 
limmass(t'/^) = limmass(i/,). By construction, we know that 



1 = id, ydivi) = ciiu -Mi)) + Zi - Aii dui. 
Since Zi — A^,- G P, we have 

j Zi-A^idi^i < j Zi -A^idvf < \\zi-A^i\\ mass{u^). 
Using Cauchy-Schwarz, and the two equations above, we get 



1 < ||v,-|| y|c,-|2 + (mass(zy,^))2. 

Now ||v,|| converges to Dist(/, O), so we find lim ||(c,-, > 1/Dist(/, O). But 
the limit of ||(c,-,f,)|| (which cannot be smaller) equals 1/Dist(/, O). Therefore, 
lim mass{v^) = lim mass(i/,) , completing the proof. □ 

To apply Theorem 5.4 to optimization problems, we will let X be the space of variations ^ 
of our given configuration and Y be the set of active constraints. Then we let f{^) 
and A^(y) be the directional derivatives of the objective function and of the constraint 
y£Y. 

In this context, a configuration satisfying condition (a) of Theorem 5.4 is called strongly 
critical, and one satisfying (b) is balanced. The theorem then says that a configuration 
is strongly critical if and only if it is balanced. 

Note that our strong criticality is indeed stronger than a simple criticality condition, 
which would say that whenever /(Q = — 1 we have Dist(A^,P) > 0, or equivalently 
that no ^ has/(0 < but A^ G P. 
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Example 5.6 With X = and Y = [0, 1] , we can set f{xi , X2) = x\ and 

A(jci , X2)(y) = 2xi \Jy-y^ + X2y 
to give an example that is critical, but not strongly critical (and thus not balanced). 

However, when F is a finite set (with the discrete topology), strong criticality is 
equivalent to criticality. For suppose whenever /(O = — 1 we have Dist(A^,P) > 0, 
but there is no uniform lower bound e > on this distance. For each y G F, we know 
that A^(y) is a linear functional on ^. Since there are only finitely many y, the graph of 
minj,gi'A^(3;) describes a polyhedron in X x R. Since the supremum over ^ G X is 
finite (we know it is nonpositive), it is achieved (at some ^ corresponding to a vertex 
of this polyhedron). But for any ^, the value is negative, so this supremum must be 
negative. 

This allows us to recover the finite-dimensional Kuhn-Tucker theorem: let X be the 
tangent space to R" at /?, let Y be the finite set of active constraints at p, and let/ 
and A be the directional derivatives of the objective function and the active constraints. 
Because Y is finite, (a) is equivalent to the definition of constrained criticality above, 
and we obtain Theorem 5.2. 

6 The balance criterion for the Gehring problem 

We now have all the tools we need to develop a balance criterion characterizing critical 
configurations for the link-ropelength problem. We start with definitions of criticahty, 
guided by our version of Kuhn-Tucker. 

Definition Suppose L is a rectifiable link with LThi(L) = r , and consider the Gehring 
problem of minimizing length subject to the constraint LThi > r. We say that L is: 

• a local minimum for link-ropelength if for all l! sufficiently C° -close to L with 
LThi(L') > r we have Len(L') > Len(L). 

• critical if for all i G VF°°(L) with b^^^iL) < we have 5+ LThi(L) < 0. 

• strongly critical if there exists some e > such that for all ^ G VF°°(L) with 
5^Len(L) = — 1, we have 6^hThi(L) < —e. 

With these definitions, we can now apply our Kuhn-Tucker theorem to the Gehring 
problem. 
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Theorem 6.1 (Balance Criterion) A link L is strongly critical for length when 
constrained by link-thickness if and only if there exists a positive Radon measure fi on 
Strut(L) such that, for every smooth vector field ^ along L, we have 

5^Len(L) = / As^d/i, 

yStrut(L) 

where As = (5Dist is the rigidity operator. 

Proof We will apply Theorem 5.4 with X := VF°°(L) and Y := Strut(L), letting 
/ := 5Len(L) be the derivative of length and A : = As be the rigidity operator We have 

||(AsO"ll = -min5pist{x,3;} 

Strut 

(when this is nonnegative). By Corollary 3.2, the right-hand side is —5'^ LThi(L), so 
that condition (a) from Theorem 5.4 is exactly strong criticality. □ 

6.1 Smoothness of critical curves 

It is unclear, a priori, how much regularity one should expect for ropelength-critical 
curves in the Gehring problem. But we can use the balance criterion to deduce 
immediately that they must have finite total curvature. 

Corollary 6.2 If a link L is strongly critical for the Gehring problem, then L is FTC . 
Proof The theorem tells us that L can be balanced: 



(5^Len(L) = / As^ d/i. 

^Strut(L) 



/Strut(L) 

But the right-hand side is a distribution of order zero on ^, since 

As^d/i < mass(/i) sup |^|. 

Strut(L) L 

Therefore, by Lemma 4. 1 , L € FTC . □ 

We can now rewrite the conclusion of our balance criterion in terms of the curvature 
force /C on L and the adjoint A^ of the rigidity operator 

Corollary 6.3 A link L is strongly critical for link-ropelength if and only if it has 
finite total curvature and there exists a positive Radon measure jj, on Strut(L) such that 

Alifi) = -/C 

as forces along L. 
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Proof The theorem guarantees that for all smooth ^, we have 

(5^Len(L) = / As^d/i. 

yStrat(L) 

By the corollary, L is FTC, so the left-hand side can be rewritten as — J^(^,d/C). 
Approximating any continuous vector field uniformly by smooth ones, we find that 

- / (e,d/C) = / AsCd/i 

JL JStrut(L) 

for all i £ VF(L), or in other words, — /C = A^ifi). □ 
We get an immediate and useful geometric corollary to this balance criterion. 

Corollary 6.4 Suppose L is critical for link-ropelength, and E C L is a subset with 
nonzero net (vector) curvature / JC{E) G M? . Then there must be at least one strut 
{e, x} with e G £ and x ^ E, and K,{E) is in the convex cone generated by the directions 
X — e of all such struts. 

Proof First note that struts from E to E contribute no net force. By the balance 
criterion, we have }C{E) = —A'^niE), and the latter is a (positive) weighted sum of 
vectors x — e. □ 

We note that this corollary is the analogue for link-ropelength of von der Mosel and 
Schuricht's "Characterization of Ideal Knots" [33, Theorem 1]. 

We next find that critical links are as well as FTC : 

Proposition 6.5 If L is strongly critical for link-ropelength, then L is . 

Proof We already know that L has finite total curvature; it is precisely when it has 
no corners, that is, when the curvature force K, has no atoms. If T± are the right and left 
tangent vectors to L at x, then /C({x}) = r+ + r_ . When K,{{x}) / 0, Corollary 6.4 
says there is at least one strut {x,y} with iy — x, /C({x})) > 0. That is, 

{y - x,T+) + {y - x,T^) > 0, 

so we must have iy — x, 7+) > or (y — x, r_) > 0. (See Figure 3.) In either case 
it follows that there exist points on L near x that are closer to y than x is, which 
contradicts the hypothesis that {x, y} was a strut. This completes the proof. □ 

The example of the tight clasp in Section 9 shows that critical links need not be 
C^'' — their curvature need not be bounded — but so far this is the worst behavior we 
can display. We conjecture that the curvature measure is always absolutely continuous 
with respect to arclength. 
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Figure 3: This curve L has a corner at x with left and right tangent vectors r_ and r+, whose 
sum is the curvature force /C({x}) there. If L is to balance, there must be a strut {x, y} with y 
in the open hemisphere (shown in light gray) of vectors with positive inner product with K.{{x}) . 
But for any y outside the normal cone (shown in dark gray), there are points near x on L that are 
closer to y than x is. Thus our {x,y} cannot be a local minimum of the self-distance function. 
This contradiction proves that a critical curve cannot have a corner 

6.2 Constraint qualification in the sense of Mangasarian-Fromovitz 

Corollary 6.3 will be the basic model for balance criteria for generalized links, and 
for links constrained by other thickness functional [2]. In some cases, including the 
link-ropelength for closed links we are treating now, we can improve on this form of 
the criterion by replacing strong criticality with criticality. This is our next goal. 

In Section 5, we defined a regular or constraint-qualified point for a finite set of 
constraints , . . . , g„ : such a point has some variation direction v such that D^gi > 
for all the active gi. By Corollary 3.2, the corresponding idea for a link L in the Gehring 
problem is the existence of a vector field ^ for which (5^LThi(L) > 0. But this is 
automatic: dilating L increases LThi to first order. 

This regularity for our problem allows us to prove that local minima are critical and that 
critical points are strongly critical. 

Proposition 6.6 A link L is critical for the link-ropelength problem if and only if it is 
strongly critical. If L is a local minimum, then L is critical. 

Proof Suppose L is a local minimum but not critical. Then for some ^ G VF°°(L) 
we have (5^Len(L) < but 5^ LThi (L) > 0. Then for small enough t > 0, the link 
L + has less length than L. This contradicts minimality unless 6^ LThi(L) = and 
thickness has decreased (but not to first order). But in this case, we can instead use the 
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rescaled deformation (LThi(L + tS)) ^ (L + t^) , for which LThi = 1 . For small t > 
these again have less length than L, contradicting minimality. 

Strong criticality always implies criticality. Conversely, suppose a closed link L 
is critical but not strongly critical. Then there exists a sequence ^, G VF(L) with 
(5g.Len(L) = — 1 and 5^ LThi(L) —>■ . Let t] be the vector field along L induced by 
dilation, scaled so that (5,,Len(L) < 1. Then we observe that 5^+g,Len(L) < for all /. 
The superlinearity of Corollary 3.3 shows that 

lim LThi > 5+ LThi > 0. 

But then for some /, we must have 

(5^+g,Len(L) < and 6+^^_ LThi(L) > 0, 

contradicting the criticality of L. □ 

Thus for closed links, a minimizer (or more generally any critical point) for the link- 
ropelength problem is strongly critical, and hence by Corollary 6.3 its curvature force is 
balanced by some strut force Ag^u. However, in our generalized ropelength problems, 
with endpoint constraints and obstacles, constraint qualification will not always hold. 
Then we will have to be careful about the distinction between criticality and strong 
criticality. 

6.3 Existence of minimizers 

We now show that each link-homotopy class contains a globally length-minimizing 
curve with LThi > 1 . 

Proposition 6.7 In a given link-homotopy class [[L]], among all curves with link- 
thickness at least 1 , there is some Lq of minimum length. 

Proof We may rescale the initial L so that LThi(L) > 2. Thus if the C" distance 
between L and a link L' is less than 1/2 then the straight-line homotopy between them 
is a link homotopy, and LThi(L') > 1 . Taking L' to be a standard smoothing of L (for 
instance, its convolution with a smooth bump function), it follows that [[L]] contains a 
C°° Unk. 

In particular, the set of rectifiable links in [[L]] with link-thickness at least 1 is nonempty. 
Let Li , L2, . . . be a sequence of such links with lengths tending to the infimal length £ 
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in this class. By the Arzela-Ascoli theorem, taking a subsequence we may assume 
that the L, converge in to a limit Lq . Since LThi is continuous with respect to the 
C° topology, and length is lower semicontinuous, it follows that LThi(Lo) > 1 and 
Len(Lo) < ^. By the remarks of the last paragraph, Lq is link homotopic to L,- for large /, 
and therefore Lq G [[L]] . Thus Len(Lo) = I and Lq is the required minimizer □ 

Since C°° links are tame, the argument above also shows the following: 

Proposition 6.8 There are no wild link-homotopy types with finitely many components. 

(This was originally observed by Milnor [25].) Thus in the work to come, we need only 
to consider tame links. 



7 Examples of critical links 

7.1 The known length-minimizing links 

In [4] , we showed that if one component of a link is linked to k others then its length 
is at least a certain constant . Although our theorem was written for the original 
ropelength problem, the proof is valid for the Gehring problem as well. Whenever a 
link can be realized with each component having length , that configuration is thus a 
length-minimizer not only when constrained by thickness but also when constrained by 
link-thickness. (These are still the only examples known to be ropelength-minimizers.) 

To any link L we can associate a graph: the vertices are the components of L, and the 
edges record which pairs are nontrivially linked. For any tree T with n edges, there is a 
unique link H(T) that is a connect sum of n Hopf links and whose associated graph 
is T. 

For many trees T with vertices of sufficiently low degree, we can realize [[H(T)]] 
explicitly with each component having exactly its minimum possible length P^ . Even 
some slightly more complicated links, like the example in Figure 4, whose graph is not 
a tree, can be realized in this way. The distance between any two linked components 
is exactly 1 . Each component in one of these minimizers is a convex plane curve 
built from circular arcs of radius 1 and straight segments. It is an outer parallel (at 
distance 1/2) to ^ shortest curve surrounding n disjoint unit-diameter disks in the plane. 
(See Figure 5.) 



Qeometry & Topology 10 (2006) 



2076 



Cantarella, Fu, Kusner, Sullivan and Wrinkle 




Figure 4: This link of six components is a global minimizer for the link-ropelength problem. 
Each component is a convex plane curve that minimizes its length given the number of other 
components it links. 




Figure 5: Here we see perimeter- minimizing enclosures of « = 1 , 2, 3 and 4 unit-diameter 
disks in the plane. The components in the known minimizing links are outer parallels to such 
curves at distance 1/2. When « = 4, the minimizer does not have a unique shape; instead 
there is a one-parameter family of minimizers. In the last shape on the lower left, there is one 
additional isolated strut, but it carries no force in the balancing measure. 

Consider the n-star T„ , the tree with a central vertex incident to all n edges. For 
n < 5 , the construction above produces a link-ropelength-critical configuration of H{Tn) 
that is known to be minimizing. We will examine the case « = 2 in detail, in light 
of our balance criterion, and then indicate how to produce link-ropelength-critical 
configurations for all n . 

Example 7.1 The link H(T2) is the simple chain of three components, shown in 
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Figure 6. In the ropelength minimizer, the two end components are circles C\ and C2, 
while the middle component is a stadium curve S. The centers of the circular arcs in S 
are points c, G C;, while the center of each C,- is a point Si ^ S. The struts are exactly 
where different components are at distance 1 . There is a strut from each point along 
each circular arc to the center of that arc (from C,- to Si and from S to c, ). There is also 
one further strut {c\ , C2} . 




Figure 6: This simple chain is known to be a minimizer for the Unk-ropelength problem, so by 
the balance criterion, its curvature force must be balanced by some measure on the struts. At the 
top, we see how the curvature forces along the circular components C, are balanced by the struts 
coming into the centers . They produce no net force on either center si . At the bottom, we see 
how the curvature forces along the semicircles of S are balanced by struts to their centers c, . 
The resulting net inward force on the c,- is balanced by an atomic measure on the one remaining 
strut {ci,C2}. 

Since we know that this configuration is length-minimizing when constrained by link- 
thickness, these struts, by Corollary 6.3 and Proposition 6.6, must support a balancing 
measure /i. Conversely, exhibiting such a measure will re -prove that this configuration 
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is critical for the link-ropelength problem, though to re -prove it is a local minimum 
would require some second-order theory. We now provide such a measure, which will 
be a useful comparison of the results of this paper against the results of [4]. 

Except for c, , each point x along the component C, is part of a unique strut {x,Si} . 
The measure assigned to struts in this "wheel" must exactly balance the curvature force 
d/C = N As{x) along C,. Because the wheel forms a complete circle, at the center 
points Si , the incoming forces from these struts cancel one another, leaving no net force. 

The situation on the stadium curve is slightly more complex. The struts from the 
semicircles of S to the points c, again balance d/C = A'^ As{x) , now for x along the 
semicircles. Unlike the previous situation, however, these measures have a resultant 
inward force of magnitude 2 at c, , directed parallel to the straight segments in the 
stadium curve. To balance these forces, the measure must have an atom of magnitude 2 
at the one remaining strut {ci , C2} . 

The measure /i we have described does balance the curvature force everywhere along 
the link, and thus demonstrates that the link is critical for link-ropelength. 

It is worth emphasizing the fact that the inner strut {c\ , C2} bears an atom of ^i. This 
stresses the point that in our Kuhn-Tucker theorem and the resulting balance criterion 
we are required to view the Lagrange multiplier ^ as a Radon measure in the dual space 
C*(Strut(L)), rather than as a density function on struts. 

Although ropelength-minimizing. Example 7.1 is not rigid, in the sense that the 
components Q can be pivoted around the points c; to be centered at any points Si on 
the semicircles of S. 

A stronger form of nonuniqueness is exhibited by the minimizing configurations [4] of 
the five-component link H{T^), with one component linked to all four others. Here the 
central component does not even have a uniquely determined shape. Instead there is 
a one-parameter family of minimizing shapes, corresponding to the deformation seen 
in Figure 5 for n = 4. Again, each of the minimizers can be balanced. (As we have 
proven, the existence of the balancing measure ji is equivalent to strong criticality for 
the ropelength problem, but it does not imply that the critical point is isolated.) 

For n > 5, we expect that similar configurations of H{T„), like the one shown in 
Figure 7 for n = 6, are again minimizing. Our balance criterion lets us show they are at 
least critical: 

Proposition 7.2 Suppose P is a convex planar n-gon with unit-length sides and 
turning angles in [0, ^'^/s] . Let Lq be the outer parallel at distance 1 from P, and 
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Figure 7: This configuration of H{T^) is critical for the link-ropelength problem. It is also 
presumably the minimizer, even though it does not minimize the length of the long component Lq 
alone. That component is the outer parallel at distance 1 to a convex planar polygon P . Each 
other component L, is a unit circle passing through a vertex of P and lying in a perpendicular 
plane. We have drawn the unit-diameter disks around these vertices, where the thick tubes 
around the L, intersect the plane inside Lq . 

let Li, . . . ,L„ be unit circles, perpendicular to the plane of P, passing through the 
vertices of P, and centered at points on Lq . Then the link L = Lq U Li U • • • U L„ is a 
configuration of H{Tn) with link-thickness 1 that is critical for link-ropelength. 

Proof As in the simple chain, each circle L, focuses a wheel of struts to its center 
point on Lq, and a measure assigning force ds to these struts balances the curvature 
force on each circle while exerting no net force on Lq . 

Let c, be the vertex of P on L, , and let 2a, be the turning angle of P there. The 
condition a, < ^^/a exactly suffices to know that no two vertices (and thus no two L, ) 
are at distance less than 1 from each other, confirming that LThi(L) = 1 . The curve Lq 
includes an arc of the unit circle around c, ; from this arc of length 2q;, a fan of struts 
converge to c, . To balance the curvature force on Lq , these struts again have measure 
equal to ds, giving a net inward force of 2 sin a,- on c,-. The remaining, isolated struts 
of L connect successive c, along the edges of P. Unit atoms of compressive force on 
these isolated struts produce exactly the outward forces 2 sin a, at c, needed to balance 
the inward forces from Lq . 
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By Corollary 6.3, the existence of this balancing measure on the struts proves that L is 
critical. □ 

For n < 5, we know these configurations for H{T„) are ropelength minimizers. For 
n > 5, the component Lq, having length n + 27r, is longer than it needs to be: at the 
expense of lengthening some other components, it could be shortened to length P„ , 
which, asymptotically, is much smaller, being 0{y/n). However, calculations we 
have done suggest that the tradeoff is not worthwhile and so the critical configuration 
described above is probably the global minimum for ropelength. 

The examples given in Proposition 7.2 — critical configurations and presumed minimizers 
for H{T„) — are quite interesting. The shape of Lq is free to move in an (n — 3)-parameter 
family; each other component is free to pivot (about its vertex of P and along one of 
the arcs of Lq), giving an additional n parameters for the shape of the whole link L. 
We also note that these examples are tight links that are not packed tightly: Consider 
the thick (unit diameter) tube around one of these configurations. As n increases, it 
occupies an ever smaller fraction of the volume of its convex hull. This should be 
compared with experiments of Millett and Rawdon [23] on this volume fraction. 

Although we have stated Proposition 7.2 above only for stars r„ , the same balancing 
works for the links H(T) based on other trees T. Each component linked to n others 
should have the shape of Lq above. We note, however, that critical links built in this 
way are not always minimizers. 

Example 7.3 Consider the tree r„ ,„ with n -\-m vertices, including one of valence n 
connected to another of valence m. The link H{Tn^,n) then has two long components, 
Lq and Li , linked to each other and to n — 1 and m — I short components, respectively. 
For large enough n and m, we can construct two thick versions of H{Tn^m), called the A 
and B configurations, as follows. 

For the A configuration, we follow Proposition 7.2 to build Lq and L\ as outer parallels 
to convex planar n- and m-gons Pq and Pi in perpendicular planes. If « = 2k, then we 
let Pq take the shape of a (^ — 1 ) x 1 rectangle capped with equilateral triangles on both 
short sides; if « = 2/c + 1 , we omit one triangle. (The precise shape is unimportant, but 
we need at least one sharp angle on each polygon.) Further, we choose the tip of such a 
triangle as the vertex of Pq corresponding to L\ (and as the vertex of P\ corresponding 
to Lq). Then it is straightforward to check that the other components stay sufficiently 
far from each other for this A configuration to indeed have LThi = 1 ; its total length 
(27r+ \){n + m). 
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However, when n and m are large enough, we can save two units of length as follows. 
Construct a tight configuration of H{Tn-\) as in Proposition 7.2, using a regular polygon. 
(Again, the precise shape is not important, but here we need a large hole in the middle 
of the polygon.) This configuration, and indeed the unit-diameter thick tubes around 
its components, is contained in a round solid torus C/„ of minor radius 3/2 and major 
radius 1 + l/(2sin ^^/ji-i). Then construct the analogous configuration of H{T,„^\) 
contained in a solid torus [/,„ . Finally, place these two pieces in space so that C/„ and Um 
form a (loose) Hopf link. (This is possible as long as the major radii are at least 3 , 
corresponding to n,m > 14.) The resulting link is the B configuration of HiTn^m)- 
Because the large Hopf link is loose, there are no struts from H{T„^\) to H{Tm-\)- 
Since each of these pieces is balanced, so is the B configuration of H{Tn,m)- 

If, as we believe, these B configurations are the ropelength minimizers, then they are 
the first ones known in which certain pairs of linked components are not in contact. (We 
note that the same must be true for the ^-component Hopf links for large n, since their 
minimum ropelength [3] is Oirv"/^). There, however, no explicit candidate minimizer is 
known. And our B configuration here has the additional property that certain linked 
pairs are not even connected by chains of touching components.) 

In all of the examples H{T) discussed above, each component is a convex plane curve 
built from straight segments and arcs of unit circles. The proven minimizers are 
minimizers in their isotopy class for the original ropelength problem [4], as well as 
minimizers in their link-homotopy class for link-ropelength. In fact, in [2] we will 
consider a family of thicknesses with varying stiffness. Each of these thicknesses is 
characterized by a stiffness A , meaning a lower bound on the diameter of curvature for 
a unit-thickness curve. Our H{T) are ropelength-minimizers for the whole family, as 
long as the stiffness A does not exceed 2 , when circular arcs of larger diameter would 
be needed. We will also develop an analog of our balance criterion for these other 
ropelength problems, and will see that all the H{T) discussed above (including those 
that are not minimizers) are critical for all formulations of ropelength where A < 2. 

7.2 Local minima for ropelength 

We do not attempt in the paper to discuss second-order behavior of ropelength near a 
critical point — in particular we have no way yet to distinguish between local minima 
and saddle points for this problem. Of course, the known minimizers must be local 
minima, and it is also easy to give critical configurations which are not local minima, as 
in Example 7.4 below. 
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Many researchers have used numerical simulations of the ordinary ropelength problem 
to look for nontrivial local minima for knots, in particular for the unknot. Such configu- 
rations have been termed Gordian unknots since they can be untangled topologically but 
not physically. Pieransky et al [29] have numerically simulated a reasonable candidate 
for a Gordian unknot, but we are very far from being able to prove its existence. 

In Example 7.3 we gave two distinct critical configurations for //(r„ ,„), and we expect 
that this will lead to the provable existence of two distinct local minima. In particular, 
our investigations lead us to predict that one cannot move from the A configuration of 
length {2iT + l)(n + m) to the suspected global minimum B without first increasing 
ropelength. This shows there must be a second local minimum; we expect, however, 
that this is not A but instead a third configuration of intermediate length. 

This connects back to Alexander Nabutovsky's original work on ropelength in higher 
dimensions and codimensions [27]. He showed using recursive function theory that, 
in those higher dimensions, a ropelength constraint often introduces new components 
into the moduli space of unknotted hyperspheres; in particular there are infinitely many 
local minima for ropelength. While for two-spheres in or for circles in there 
are presumably no such minima, we do expect there must be infinitely many Gordian 
unknots in M? . Our two critical configurations of HiTn^m) are perhaps a first step toward 
proving this. 

7.3 Elastic tension energies 

All of the links presented above are critical or minimizing for the sum of the lengths 
of their components. This is a beautiful functional, but it is physically somewhat 
unrealistic: elastic ropes should minimize a quadratic functional of the form 

y^fl!, (Len(LO - 

i 

where a, > is the elasticity and li the rest length of the /-th component. Criticality 
for this functional is equivalent to that for ^ ti Len(L,) where the tension ti is tj = 
2a,(Len(L,) — li). Assuming these tensions are nonnegative (that is, that no component's 
length is less than its rest length) our balance criterion extends immediately to handle 
this case: the strut force A^{ix) must balance the tension-weighted curvature force 

In the known minimizing links, such as the simple chain, each component separately 
achieves its minimal possible length. Thus these examples also minimize all elastic 
energies with nonnegative tensions > 0. 
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This behavior, however, seems rather exceptional. The examples in Proposition 7.2 do 
not minimize all such functional. In particular, if the tension in the long component is 
large enough, it will shrink to length while some of the shorter components gain 

length. 

Also in the Borromean rings, if the three components have different tensions, the 
configuration we describe below (Section 10) would no longer be critical. Similarly, 
clasps (Section 9) in which the two ropes have different tensions again have new critical 
configurations. In [37] we describe in detail the shapes of these asymmetric clasp 
curves, as well as their appearance in more complicated clasp-like links even when 
tensions are equal. (Note that link-thickness was called Gehring thickness there, as in 
early drafts of this paper.) 

7.4 Nonembedded critical links 

To illustrate the differences between the Gehring problem and the original ropelength 
problem, we now give some examples of a different flavor: critical configurations that 
are nonembedded and thus have infinite ropelength in the original sense. 

Any knot is of course link-homotopic to the unknot. The link-ropelength minimizer 
degenerates to a point (of length zero). The same happens for any component of an 
arbitrary link that is link-homotopically split from the rest of the link. 

Milnor showed that, up to link homotopy, links of two components are classified by 
their linking number [25]. 

Example 7.4 When the linking number is zero, the components split, and the link- 
ropelength minimizer degenerates to have length zero. We can, however, also describe 
another critical configuration for this unlink: one component degenerates to a point p 
while the second is a unit circle centered at p. This is clearly an unstable critical point: 
obvious deformations can decrease the ropelength to second order. 

The case of linking number 1 is close to Gehring's original problem: the minimizer is 
the same Hopf link built from round circles. (This case fits in the class H{T„) considered 
above.) For larger linking number, we can use Corollary 6.3 to exhibit many critical 
configurations as follows: 

Example 7.5 For linking number mn there is a critical configuration L„, „ consisting 
of the minimizing Hopf link with one component covered m times and the other covered 
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n times. Its total length is thus 27r(m + n). There are other critical configurations, 
sometimes shorter. For example, each component can be a figure-eight built from two 
tangent circles. Figure 8 shows a configuration like this with total length 27r(m + n) and 
linking number mn — m\n\. The best configurations we know for linking number 17, 
for instance, use {m, n) = (6, 3) or (4, 5). Assuming configurations like these are the 
minimizers for two-component links, they give examples where the set of minimizers is 
disconnected (since we can interchange the two components, or reorder the way one 
component covers its figure-eight). 



Figure 8: In this configuration of two curves from Example 7.5, each circle is covered m,- or n, 
times, as labeled. If m = mi + m2 and n — ni + n2, then the curve has total length 27r(m + «), 
and linking number mn — mini. It is constrained-critical, though often not minimal, for the 
Gehring problem in the link-homotopy class defined by its linking number. 

None of these configurations is embedded, so they are not critical points for the original 
ropelength problem: as expected, the extra freedom in the Gehring problem sometimes 
allows for shorter solutions. As a further example, consider the (2, 4)-torus link, with 
linking number 2 . We have computed the presumed ropelength-minimizer numerically, 
as in [35]. The results are shown in Figure 9; this solution is longer than the covered 
Hopf link L2,i (the presumed link-ropelength-minimizer) and is not even critical for the 
Gehring problem. 

For links of more than two components, linking numbers do not suffice to distinguish 
link-homotopy types; we must also consider Milnor's /i -invariants [25, 16]. For 
instance, the Borromean rings, with no nonzero linking numbers, belong to a nontrivial 
link-homotopy class because they have /i -invariant equal to 1 . 

Numerical experiments performed with Brakke's Evolver (compare [35]) suggest that 
the minimizing Borromean rings for the link-ropelength problem should consist of three 
congruent curves in perpendicular planes. In [4], we described such a configuration built 
from circular arcs of radius 1 . Unfortunately, Corollary 6.3 shows this configuration 
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Figure 9: This picture shows a numerically computed minimizer for the original ropelength 
problem on the (2,4)-torus link. Because it has a strut between two points on the same 
component (shown center, where the darker tube contacts itself) that carries nonzero force, it is 
not balanced for the Gehring problem considered here. It is longer than L2.1 , the Hopf link with 
one component doubly covered, which we conjecture is the minimizer for link-ropelength in 
this link-homotopy type. Notice that both of these configurations break the symmetry between 
the components of the link, so we expect to also find a (longer) critical configuration where the 
two components are congruent. 

is not even critical for length when constrained by Unk-thickness. In Section 10, the 
culmination of our paper, we will explicitly describe a very similar configuration of the 
Borromean rings, which we prove is critical and believe is the minimizer. 

However, in order to solve for these Borromean rings, we must first consider a simpler 
interaction between two ropes: the clasp that occurs when one rope is pulled over 
another. Describing this will require a notion of generalized links. 

8 Generalized link classes 

Although some of our definitions have applied to arbitrary curves, so far we have been 
treating only ordinary (closed) links. We now want to consider generalized problems 
involving curves with endpoints. To get meaningful link classes in this setting, we must 
include constraints for the endpoints and obstacles for the link. 

Definition A generalized link L is a curve L (with disjoint components L\, . . . ,Lfj) 
together with obstacles and endpoint constraints. In particular, each endpoint x G 9L is 
constrained to stay on some affine subspace Mx C , which can have dimension , 1 
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or 2 . Furthermore, there is a finite collection of obstacles for each component L,- of the 
link. Each obstacle 

{/p G : gijip) < 0} 

is given by a function gjj with as a regular value. By calling them obstacles, we 
mean that L, is constrained to stay in the region where minyg,y > 0. 

While we could allow even more general endpoint and obstacle constraints, this version 
fits nicely with our overall setup, and allows for all the specific examples we have in 
mind. 

Definition Suppose L = [j^ L, is a generalized link, with obstacles gij and endpoint 
constraints Mv. Then its link-homotopy class [[L]] is the set of all links L' that are 
link-homotopic to L through links where each component avoids its obstacles and 
maintains its endpoint constraints. (As before, in a link homotopy, each component 
of L can intersect itself but not the others.) 

This definition is comparable to our previous definition for closed links (Section 2); 
as in the discussion at the end of Section 6, we may restrict our attention to tame link 
classes. 

Given a generalized link L, only variations preserving the endpoint constraints should 
be allowed. A vector field ^ € VF(L) is said to be compatible with the constraints 
if it is tangent to at each endpoint x G 5L. We write VFc(L) for the space of all 
compatible vector fields. 

Given a set of obstacles gij < and a link L = |J L,- , we write 

0{L) := minmin p,7(x). 

i.j xdLj 

Then L avoids the obstacles gij if and only if 0{L) > . We define the wall struts of L 
by 

Wally(L) := Li n {gij = 0}, Wall(L) := |J Wally(L). 

'J 

This incorporates those parts of L on the boundary of the obstacle, but is not strictly 
speaking a subset of L since one point x € L; might be in several of the Wall,y . When 
0(L) = 0, by Clarke's Theorem 3.1 we have 
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Again, we collect the various derivatives appearing on the right-hand side into a rigidity 
operator Aw : VFc(L) — > C(Wall(L)) on wall struts, given by 

AweW := (6, Vgij) 
when X E Wally . Its adjoint is then 

(8-1) [^dA*^{fi)=[ Awed/x = V/ {i,,Vgij)dfi(x). 

JL JWall(L) •• J.veWall,j(L) 

'5/ 

We also have corresponding definitions for locally minimal, strongly critical, and critical 
configurations of L : 

Definition We say that a generalized link L is a local minimum for length when 
constrained by LThi if we have Len(L') > Len(L), for all sufficiently C"-close links L' 
with the same obstacle and endpoint constraints and with LThi(L') > LThi(L). We say L 
is strongly critical (respectively, is critical) for minimizing length when constrained by 
LThi if there is e > such that for aU compatible smooth ^ with 5^ Len = — 1 , the 
quantity 

min (5+LThi(L), (5+0(L)) 
is at most — e (respectively, is negative). 

As in our discussion of Kuhn-Tucker at the beginning of Section 5, these notions 
will be equivalent only under a regularity assumption corresponding to the constraint 
qualification of Mangasarian and Fromovitz [20] : 

Definition A generalized link L is YXhi-regular if there is a thickening field, meaning 
a smooth compatible rj for which 5^ LThi(L) > and 5^0{L) > 0. 

Note that, while we require rj to strictly increase LThi and to move away from the 
obstacles, both to first order, there is no corresponding requirement for the endpoint 
constraints, since they are linear equality constraints instead of nonlinear inequality 
constraints. 

We can now prove a generalization of Proposition 6.6: 

Proposition 8.1 If a generalized link L is a LThi -regu7ar local minimum when 
constrained by LThi, then L is critical. Also, if L is LThi -reguJar and critical when 
constrained by LThi, then it is strongly critical. 
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Proof The regularity of L means there exists a thickening field rj € VFc(L) . We may 
assume (5,;Len(L) > for otherwise L is neither minimal nor critical; we then scale rj 
so that 5^Len(L) < 1 . 

Suppose that L is a local minimum but not critical. Then for some compatible vector 
field ^ we have 5^l^n{L) < while (5+LThi(L) > and 5^0(1) > 0. For small 
f > 0, consider the links Lt = L + t{£, + er/). Then 
dLen(L,) 



^ = 5^Len(L) + e(5^Len(L). 

We choose < e < — (5^Len(L)/5^Len(L), so this derivative is negative at time 0. Thus 
for small t, the Lf have length less than Len(L), contradicting minimality if they obey 
our constraints. But 



dLThi(Lf) 



d^^ 



> 0. 



>0, 

f=0 



t=Q ' dt+ 

and the endpoint constraints are linear, so the links L, meet all our constraints for small 
t > 0. 

Now suppose that L is critical without being strongly critical. Then there exists a 
sequence of compatible vector fields S VFc(L) with (5g.Len(L) = — 1 but with either 
(5^LThi(L) — > or 6^,0{L) —>■ 0. Then we observe that (5^_|_^;Len(L) < for all /, 
while by Corollary 3.3 either 

lim6\, LThi > 6t LThi > 
or lim(5+ O > 5+0 > 0. 

Taking / large enough that one of these quantities is positive, we get a contradiction to 
the criticality of L. □ 

So far, this development has paralleled that of Section 6; we now diverge from our 
previous course. Earlier, we saw that every closed link is LThi-regular: rescaling 
always provides a thickening field. In the generalized setting, this is no longer the case. 
Thus minimality no longer implies criticality. 

Example 8.2 To give a specific example, rotate the constraints of Example 5.1 around 
the z-axis to give obstacles gi = (x^ + — l)^ — z < and g2 = z < for an 
unknot L. The unit circle in the xy-plane is on the boundary of both obstacles, and 
is clearly the minimum-length configuration in its homotopy class. However, it is not 
critical: shrinking it toward the origin will reduce its length to first order; the constraint 
gi > is now violated, but not to first order. 
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Further, criticality and strong criticality may be different: if we allowed infinitely many 
obstacles, we could construct critical, but not strongly critical links by following the lead 
of Example 5.6. (If we do not allow infinitely many obstacles, then an open question 
remains: is strong criticality a stronger assertion than criticality?) 

Example 8.3 To justify our emphasis on strong criticality (rather than restricting our 
attention to regular, critical links) we also note that it is easy to construct strongly 
critical links that are not regular; simply take L to be the unit circle in the xy-plane, 
with constraints gi(x, y,z) = + — i (so the excluded region is the infinite cylinder 
around the z-axis) and g2 = —g\ ■ This link is trapped on the cylinder = = ^2. so 
it has no thickening field. On the other hand, it is clearly strongly critical. 

Now we are ready to extend our balance theorem to the generalized setting. We will 
accommodate the endpoint constraints by restricting our attention to compatible vector 
fields. Our other constraints are then Dist > 1 on L^^^ and gjj > along L, . The set Y 
of active constraints then consists of the struts together with the wall struts. 

Theorem 8.4 A generalized link L is strongly critical for link-ropelength if and only 
if there is a positive Radon measure /x on Strut(L) U Wall(L), such that 

-K = (As eAw)> 

as linear functionals on VFc(L). This means that —K and (As © Aw)*^ agree as forces 
along L except at endpoints x ^ dL, where they may differ by an atomic force in a 
direction normal to Mx . 

Proof This is again a straightforward application of our Theorem 5.4, using 

X = VFc(L), F = Strut(L) U Wall(L), f = 6Len, A=AseA^f^. □ 

Remark 8.5 Remember that /C has been defined to include an inward-pointing atom 
at each endpoint x G dL. We can ignore these, however, when applying this theorem, 
as long as the link L meets each endpoint constraint normally. We know of no 
examples of critical links where this is not the case. 

The regularity described in Corollary 6.2 and Proposition 6.5 carries over to generalized 
links: 

Proposition 8.6 If the generalized link L is strongly critical for link-ropelength, then L 
is FTC and . 
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Proof The proof follows that of Corollary 6.2 and Proposition 6.5. From equation (8-1 ), 
we find that has distributional order zero just like A^/x, so L S FTC follows 
immediately from the balance criterion of Theorem 8.4. 

Now suppose L is not but instead has some comer x with lC{x} 7^ 0. By 
Theorem 8.4, this curvature force is balanced by struts and wall struts. So there is at 
least one strut or wall strut acting on x in a direction with negative inner product with 
/C{x}. In the case of a strut {x, y}, we refer again to Figure 3: some points near x 
along L would be nearer to the endpoint y. But similarly, in the case of a wall strut, 
we have (/C, Vg,y) < 0, but this means that some points near x along L violate this 
obstacle constraint. In either case, we get the desired contradiction. □ 

To understand the interplay between struts and wall struts, we now offer a simple 
example of a generahzed link L with nonempty boundary which is balanced, needing 
nonzero force on the wall struts. 

Example 8.7 Cut the simple chain of Figure 6 by parallel planes through s\ and S2 
with normal vector c\ — C2, and let L be the part of the chain lying between the two 
planes. This generalized link includes two semicircles with endpoints normal to the 
planes, and also the inner stadium curve, which is tangent to the planes at and S2- 
We let the planes bound an obstacle, forcing L to stay between the planes, and we use 
them also as endpoint constraints. Then L is balanced: though the semicircles now 
exert a net outward force on s\ and S2, this is balanced by wall struts at these points. 
And the internal balance for the stadium curve remains the same. 

9 The tight clasp 

The tight configurations of Section 7 were the simplest closed links we could imagine: 
the Hopf link, and various connect sums of Hopf links in which each component is still 
a convex plane curve. But there is an even simpler interaction between two ropes, the 
clasp formed when one rope is pulled taut over another, as at the junctions of a woven 
net, or when a bucket is lifted from a well by passing a rope through its rope handle. 
We can model a single clasp as a generalized link with endpoint constraints. 

To define the simple clasp, fix two parallel planes P and P at least 2 units apart. Then 
take two unknotted arcs 7 and 7 that lie between the planes, with the endpoints of 7 
constrained to lie in P and those of 7 in P. Let the half space bounded by P that does 
not include P be an obstacle for the component 7 , and vice versa, and select the isotopy 
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class of such links shown in Figure 10. This is the class where closing each arc in the 
plane of its endpoints would produce a Hopf link. 



Figure 10; The simple clasp has two components, one attached at both ends to the ceiling and 
the other to the floor, linked with one another as shown. The configuration shown, with each 
component consisting of a semicircle joined to two straight segments, is neither critical nor 
minimal. 

It is natural to assume that the minimizing configuration for this problem would consist 
of semicircular arcs passing through each others' centers, together with straight segments 
joining the semicircles to the constraint planes, much like the Hopf chain of Example 7.1. 
But this naive clasp is not balanced: each semicircle focuses its curvature force on 
the tip of the other, and there is no way to balance these forces (as the isolated strut 
carrying an atom of compressive force did in the Hopf chain). The naive clasp is thus not 
minimizing, though we will see it is very close: the critical configuration we construct 
here is only half a percent shorter. 

Example 9.1 Suppose the horizontal planes P and P in the definition of the simple 
clasp are taken instead to be only one unit apart. Consider the configuration where the 
curves 7 and 7 are semicircles in perpendicular vertical planes. The curvature of each 
semicircle can be balanced by uniform strut tension, transmitting a net vertical force to 
the tip of the other semicircle. That vertical force can be balanced by a wall strut at 
each tip. Therefore, this configuration is critical for link-ropelength. 

In the case of interest, where P and P are far apart and there are no wall struts to 
balance the tips, we must look harder for a solution. We will now construct critical 
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configurations, constrained by the link-thickness LThi, for the simple clasp problem 
and for a family of related problems where the ends of the ropes are pulled outward as 
in Figure 11. These solutions minimize length under natural symmetry assumptions, 
and we believe they are the global minimizers even without imposed symmetry. Below 
in Section 10, we will construct a critical configuration of the Borromean rings that 
contains portions of these clasp curves. Thus, a thorough understanding of these 
generalized links will aid us in understanding that more complicated closed link. 

9.1 Symmetry conditions and a convenient parametrization 

We describe configurations of the clasp where the two components are congruent plane 
curves, lying in planes perpendicular to each other and to the constraint planes. To 
fix these symmetries in coordinates, let the constraints be the planes z = ±C, and 
let the component 7 lie in the xz -plane while 7 lies in the -plane. The clasp has 
mirror symmetry across each of these planes (preserving each component). It also has 
a symmetry interchanging the two components, which we denote p ^ p, given by 
fourfold rotation about the z-axis together with reflection across the xy-plane. These 
symmetries generate a point group of order eight in 0(3) whose Conway-Thurston 
orbifold notation [9, 10] is 2 * 2. Algebraically it is isomorphic to D4. 

The argument we present below to derive the critical clasps for the Gehring problem can 
easily be extended to show these are the unique critical configurations among curves 
with this 2*2 symmetry. We omit the details, however, because we know of no way to 
show that the overall minimizers must have this symmetry. If one could prove this, it 
would then follow that our clasps are the minimizers. 

Our symmetry assumptions mean that the clasp is described by the shape of half of the 
component 7, from its tip along the z-axis into the x > half-plane and up to the 
plane P. This consists of a curved arc near the tip joined to a straight segment near P. 
Since the curved arc is strictly convex, we can parametrize it by the angle Lp made by its 
tangent vector above the horizontal, as in Figure 1 2. In fact, we will use the sine of this 
angle, m = sin as our parameter. Thus in the simple clasp, for u G [0, 1] we write 

7(±m) = (±x(m),0,z(m)), 
7(±m) = (0, ±x(m), -z(m)). 

Elementary calculations show the following: 
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Lemma 9.2 For a convex curve 7 in the xz-plane, parameterized by the sine u of its 
direction G [— |, |], the arclength s satisfies 

ds = sec (fdx = esc ipdz = — 1^=^= , 

KV 1 — 

where the curvature k is given by 

dip du 
ds dx 

For the simple clasp described above, each component turns a total of 180° , meaning 
that u ranges from — 1 , through at the tip, to 1 . We can also consider more general 
clasp problems where the four ends of rope are not vertical (being attached to horizontal 
planes) but instead are pulled out at some angle (being attached to tilted planes). 

Given < r < 1 , we define the r -clasp to be a problem like the simple clasp where 
the arc 7 starts at m = — r and then turns through angle 2 arcsin r to reach u = t. Our 
critical r-clasps have the same 2*2 symmetry as the simple clasp. To put the r-clasp 
into our framework of generalized links, we constrain the four endpoints to four planes, 
each making angle arcsin r with the vertical, as in Figure 1 1 . The complement of the 
wedge formed by the planes containing the endpoints of each arc acts as an obstacle for 
the other arc. The simple clasp is the r-clasp with r = 1 , where the wedges degenerate 
to halfspaces. 



9.2 Struts between perpendicular planes 

Whenever two curves in perpendicular planes are connected by a strut, elementary 
trigonometry gives us first order information about the curves at both endpoints. We 
state a general lemma, which we will use here for the clasp and again for the Borromean 
rings. 

Let Pi and P2 be two planes meeting perpendicularly along a line £, and let 7,- C Pj 
be two components of a link. At a point p, G 7, , we write x, for the distance from p, 
to £, and m, for the cosine of the angle between £ and the line tangent to 7, at p,-. These 
quantities generalize the x and u of Lemma 9.2 above. 

Lemma 9.3 Let 71 and 72 be two components of a link L, lying in perpendicular 
planes. Suppose there is a strut {p\,p2] of length 1 connecting these components. 
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arcsin r 




Figure 11: In this variant of the simple clasp problem, the endpoints of the two ropes are 
constrained to lie in four planes whose normals make angle arcsin t with the horizontal. The 
parameter m = sin ranges from — r to t along each arc, as shown at the end of the top right 
arc. If extended, the four planes shown would form the sides of a tetrahedron. Each arc is 
constrained to lie in the wedge formed by the planes containing the endpoints of the other arc. 



Then in the notation of the previous paragraph we have < x, < m, < 1 , and any two 
of the numbers xi , X2, mi , M2 determine the other two, according to the formulas 

Jl ,,2ri ,,2\ 



u} 



2 ^ X/ ^ Ufjl - Up 

' 2 1 2 2 

Uj 1 — UfUj 



2 1 2 ' 



where j ^ i. 



Proof Picking cartesian coordinates such that i is the z-axis and P, are coordinate 
planes, we find the strut difference vector pi — P2 is (xi,X2, Az), for some number Az. 
Since this strut has length 1 and is perpendicular to each 7, , we have 

Az2+x?+xf = l, Az 



Simple algebraic manipulations, eliminating Az, lead to the equations given. □ 

Note that the condition x, < a, is exactly the condition that the unit normal circle 
around p, intersects Py; the two points of intersection are mirror images (across Pi), 
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with the same xj and uj values. Also note that we don't need to have 7; C Pi in the 
lemma; it suffices that 7, be tangent to P, at . 

Whenever we have a pair of curves in perpendicular planes, which stay a constant 
distance 1 apart, we can apply this lemma everywhere along the curves. Each curve 7; 
is determined as the intersection of the plane with the unit-radius tube around the 
other curve 7^. This will be the situation for the clasp. 

9.3 The balancing equations for the clasp 

By Theorem 8.4, in a critical clasp the curvature force along 7 must be balanced by 
struts to 7. In particular, almost every point (indeed, since the set of struts is closed, 
every point) 7(m) along the curved arc of 7 must have a strut to some point 7(m*). 
Then by symmetry we actually have what we call 2-to-2 contact: there are struts 
from 7(±m) to 7(ibM*). Here the two points 7(ibM*) must be the intersection of the 
unit normal circle around 7(m) with the -plane, implying that u* G [0, 1] is uniquely 
determined for each u. We will refer to 7(11) and 7(a*) as conjugate points on the 
r-clasp. Lemma 9.3 applies to any pair of conjugate points, with u\ = u, U2 = u* and 

Xi = x{Ui). 

Lemma 9.4 Suppose 7 is a plane curve, symmetric across a line £ in the plane. 
Consider the net curvature force of a mirror image pair of infinitesimal arcs of 7 . This 
acts in the direction of the line £, v^^ith magnitude 2\du\. Here the function u is defined 
along 7 as the cosine of the angle ip between I and the tangent line to 7. 

Proof One infinitesimal arc has net curvature force kN As = N dip. When this is added 
to the mirror image force, only the component along i survives. We get magnitude 

2| sinV'd'i/'l = 2|dM|. □ 

Suppose now we have a symmetric configuration of the clasp where the curved arcs of 
the two components stay a constant distance 1 apart. By symmetry we get the 2-to-2 
strut pattern described above. Assuming the straight ends of each component meet the 
constraint planes perpendicularly, our balance criterion Theorem 8.4 says that strong 
criticality is equivalent to the statement that the net vertical curvature force exerted 
by the arcs at 7(±a) balances that of the conjugate arcs at j(±u*). That is, using 
Lemma 9.4, for a critical clasp we must have |dM| = \du* \ , meaning that either u — u* 
ox u + u* is constant. 
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If M — M* were constant, by symmetry it would be zero, and our equations would describe 
a pair of half-ellipses, with horizontal major axis \/2 and vertical minor axis 1 . On 
these curves, corresponding points 7(m) and 7(m) are always at distance 1 from each 
other, but these are maxima for the distance between components, rather than minima. 
This configuration has LThi < 1 , and is not LThi-critical: the pairs {7(m), 7(m)} are 
not struts. 

Instead we must have that m + m* is constant. To find the constant, note that on the 
r-clasp, the tip of 7 (at m = 0) is joined by a strut to the end of 7 (at u* = r); thus 
u + u* = T. This equation holds when < m, m* < r; to allow for negative values 
(parametrizing the whole clasp curve) we write 

\u\ + \u*\ = T. 

We can now give an explicit description of our critical r-clasp: 

Theorem 9.5 Let r G [0, 1], and let 7 = 7,- be the curve in the xz-plane given 
parametrically for u € [— r, r] by 

ly/l-{T-\u\)^ 



X = Xr{u) 



yJ\-u\T -\u\f'' 



, dz , f u du 

z = Zt{u) := / —dx 



dx J - l^riu) ' 

;i_„2(^_ 1^1)2)3(1 1^1)2) 

where nJu) :=- -—^ , , , , ^ — 

1 -(r- |M|)2 + (r- |m|)|m|(1 -m2) 

and the constant of integration for z is chosen so that 

z(0)+z(t) = -\/l -r2. 

Then the union of 7 with its image 7 under the symmetry group 2*2 described above 
is a T -clasp that is strongly critical for link-ropelength. The curvature of 7 is Kr{u) 
above, and the total length of the curved part of 7 is 

du 



T k^(m)V1 — 

Proof The proposition follows from the foregoing discussion, after substituting 
u* = T — \u\ into the equations of Lemma 9.3, and using Lemma 9.2. To get the 
constant of integration for z, we note that the strut from 7(6) to 7(r) has height given 
(as in the proof of Lemma 9.3) by 

Az = Vl -XriOf-Xrirf = Vl - 0-r2. □ 
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Although the formulas we have given for Zt{u) and for arclength both involve hyperel- 
liptic integrals not expressible in closed form, it is straightforward to integrate them 
numerically; we have plotted our critical configuration of the simple (r = 1 ) clasp in 
Figure 12. 



Z 




T 

Figure 12: This is an accurate plot of the critical simple clasp 7 given by Theorem 9.5. Here 
u — smip ranges from — 1 to 1 over the curved portion of 7 . The tip 7(0) of the other 
component is shown above 7 on the z-axis, along with the (dotted) circular cross-section of the 
tube of unit diameter around 7. The curved dotted lines extending down from the sides of this 
cross-section are the lines of contact between the shaded tube around 7 and the front half of the 
tube around 7. Symmetric Hnes of contact extend behind the shaded tube, reahzing the 2-to-2 
contact pattern we have described. Finally, we see a small gap between the tubes, explored in 
more detail in Figure 14. 



As we mentioned in the introduction, Starostin has given [34] an independent derivation 
(using a form of balancing for smooth curves) of these same r-clasp configurations (as 
well as the family of stiff clasps we will consider in [2]). Starostin does not prove that 
these configurations are critical for link-ropelength. 
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9.4 The geometry of the tight clasp 

We now examine the curvature and other geometric features of the critical clasps for 
the Gehring problem that were given in Theorem 9.5. Each component of the critical 
r-clasp is a C' join of four analytic pieces: a straight segment, then 7[— r, 0], then 
7[0, r] , and finally another straight segment. Where the curved arcs join the straight 
segments at m = itr , the curvature k{u) approaches 1 ; at these points, our critical clasp 
agrees to second order with the naively expected circular arcs. 

The maximum curvature k;(0) = 1 / Vl — occurs at the tip. For r < 1 , this is finite, 
and our r-clasp is C''^ . But for r = 1 , the curvature blows up (like at the tip. 

In Figure 13 we plot the curvature k{u) for this simple clasp. The curve is C^'^/^ (and 
is also in the Sobolev space VK^' for all e > 0) but has no higher regularity. 

3 




0.5 1 1.5 

s 

Figure 13: The graph shows the curvature k of the tight 1 -clasp as a function of arclength. The 
curvature blows up at the tip: this curve is only C' "^^ . The unit-diameter thick tube around 
the curve forms a cusp near the tip, when the curvature exceeds 2 . From the tip, curvature 
decreases rapidly to its minimum, and then increases again to the limiting value of k = 1 at the 
end. Thus the clasp curve, at its end, agrees to second order with the naively expected unit circle 
around the tip of the other component, as is suggested in Figure 12. (For t = 1 , as illustrated 
here, the curves agree even to third order.) 

In Proposition 6.5, we proved that critical curves for link-ropelength are . It would 
be interesting to find out whether all such critical curves are C^''/^ ; perhaps the simple 
clasp exhibits the worst possible behavior. 

In Example 7.5, we saw critical curves constrained by link-thickness which fail to have 
positive thickness in the ordinary sense of [4] because one component is nonembedded. 
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Figure 14: We see three views of tlie gap ciiamber between the two tubes in the tight clasp with 
T — 1. On the left, we see an exploded view with the two tubes and the gap chamber floating 
between them. In the medium closeup in the center, we see the chamber in place between 
the (now transparent) tubes. On the right, we see an extreme closeup of the center of the gap 
chamber. Its height at the center (about 0.05639) is the distance between the tubes at the tips of 
the clasp. The grid in the center and right pictures is a square grid projected from the xy-plane. 
On the right, we see a tiny ridge running from left to right along the surface of the chamber; this 
is a cusp formed by the folding of the tube surface that happens when the curvature of the clasp 
rises above 2 (compare Figure 13). We do not know whether this gap chamber forms in clasps 
of physical rope; it would be very interesting to find out. 

The simple clasp fails to have positive thickness for a different reason: its curvature 
is unbounded. In [2] we will consider a family of thickness measures with a variable 
stiffness parameter A . In these measures, a unit-thickness curve has curvature bounded 
above by 2/a. For any nonzero A, it follows that the critical simple clasp must be 
different from the tight clasp here for the Gehring problem, and must instead include an 
arc of this maximum allowed curvature. 

One of the most interesting features of the clasp is the gap between the two components 
of the clasp. The distance between the tips of 7 and 7 is z(t) — z(0) + \/l — (written 
in this way to be independent of the constant of integration for z). This is an increasing 
function of r, close to 1 when r is small, but increasing to 1.05639 at r = 1 . Thus, in 
the simple clasp, the gap between the thick tubes around the two components at their 
tips is almost 6% of their diameter. 

These thick tubes contact each other at the midpoints of the struts. Topologically, the set 
of struts forms a loop. Their midpoints form a loop in space with four vertical cusps — the 
line of contact of the two tubes — as seen in Figure 12 and Figure 14. Alternatively, we 
can plot the loop of struts as pairs of arclength coordinates on the two components, as 
in Figure 15. The solid tubes divide the rest of the ambient space into two regions: one 
infinite component around the outside of the clasp, and one small chamber sitting in the 
gap between the tips, shown in Figure 14. To give a sense of scale, the gap chamber 
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Figure 15: The graph shows the strut set for the tight 1 -clasp, where each strut is plotted 
according to the arclength of its ends on the two components of the clasp (measured from the 
tip J = to the shoulders at s = i ^ 1.58944). There is a closed loop of struts, with four cusps 
at the tips and shoulders of the clasp arcs. We hope this explicit strut set will help in verifying 
the accuracy of numerically computed strut sets for ropelength minimizers, such as those in [5]. 

has a substantial surface area of about 1.10, equal to the area of a section of tube of 
length more than i/3- However, the chamber is very thin, resulting in a volume of 
only 0.01425. 

9.5 Length comparison with the naive clasp 

Earlier, we described the naive circular configuration for the simple clasp. Similarly, in 
what we call the naive r-clasp, each component is built from straight segments (normal 
to the constraint planes) and a unit-radius arc (of angle 2 arcsin r and centered at the 
tip of the other component). As we saw for r = 1 , this configuration is not critical: 
there is no way to balance the forces concentrated on the tips, unlike in Example 7. 1 
and Example 9.1, which had extra struts. 

Our critical r-clasps (which we expect are the global minima for length) are indeed 
slightly shorter than the naive configurations. The total length of a clasp depends, of 
course, on the position of the bounding planes. Thus to compare the lengths of the 
naive clasp and our critical clasp in a meaningful way, we introduce the notion of excess 
length. The infimal possible length of a r-clasp with no thickness constraint is easily 
seen to be four times the inradius of the bounding tetrahedron. (In the case r = 1 this 
is twice the thickness of the bounding slab.) The excess length of any given clasp is the 
amount by which its length exceeds this value. 
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For r = 1 , the naive clasp has excess length 2tt — 2, since two unit semicircles replace 
two straight segments of unit length. Numerical integration reveals the excess length 
of our critical 1 -clasp to be 4.262897 (accurate to the number of digits shown). It is 
thus about 0.020288, or almost half a percent, shorter. In general, the excess length of 
the naive r-clasp is 4arcsinT — 2r, while the excess length of our critical r-clasp 
equals the total length of the curved parts minus 2r times the inter-tip distance. The 
maximum percentage savings, about 0.518%, occurs for r k, sin(80°). 



10 The Borromean rings 

The original Gehring link problem was solved by the Hopf link made from a pair 
of circles through each other's centers. We have already generalized this to a three 
component link in one way: the simple chain made from circles and stadium curves of 
Section 7. But the simple chain is just a connect sum of Hopf links, and so the minimizing 
configuration shares much of its geometry with the original Gehring solution. 

We now construct a proposed minimizer for a more interesting Gehring problem — the 
BoiTomean rings (see Figure 16). Among the three prime six-crossing links of three 



Figure 16: Our critical configuration Bq of the Borromean rings, shown with thick tubes of 
diameter 1 . This configuration is very slightly shorter than the piecewise-circular version in [4]. 
As in that version, the core curve of this tube has discontinuous curvature, for instance at the 
"jump point" where the curve switches from convex to concave. 

components, the Borromean rings form the one which is Brunnian, meaning that if 
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any one component is removed the remaining components are unlinked. Milnor's /i- 
invariant classifies three-component Brunnian link-homotopy types, and the Borromean 
rings are the first nontrivial example. 

In this section, we describe (Theorem 10.2) a critical configuration Bq of the Borromean 
rings, shown in Figure 16 and Figure 17. Numerical simulations with Brakke's 



Figure 17: Two further renderings of the critical configuration Bq for the Borromean rings 
reveal more of the structure. The image on the left, showing thin tubes of diameter 0.315, 
viewed along an axis of threefold symmetry, has been adopted as the logo of the International 
Mathematical Union. On the right, in a still from the video [15], we see even thinner tubes 
inside transparent thick tubes. 

Evolver [1] suggest that this configuration Bq is in fact the ropelength minimizer for the 
Borromean rings. We will see below that the curvature of Bo stays below 1.534; this 
means (as we show in [2]) that Bq is also a critical point for length when constrained by 
the ordinary thickness measure of [4] instead of by link-thickness. In [4], we described a 
similar configuration 62 of the Borromean rings, built entirely from arcs of unit circles. 
Theorem 6.1 shows that B2 is not critical, and we compute that Bq is 0.08% shorter. 

10.1 Symmetry and convexity 

Our configurations Bq and B2 of the Borromean rings are quite similar, and in particular 
have the same symmetry and convexity properties, which we now define. The three 
congruent components lie (respectively) in the three coordinate planes; reflection across 
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any one of these planes is a symmetry of the hnk preserving each component. A further 
symmetry, which cycUcally permutes the three components, is given by 120° rotation 
about the (1,1,1) axis; we write this rotation as 

p ^ p ^ p p. 

These symmetries generate the pyritohedral point group of order 24 in 0(3) whose 
Conway-Thurston orbifold notation [9, 10] is 3*2. Algebraically it is isomorphic 
to ib44 , and in cartesian coordinates it is most naturally seen as the wreath product 

{±1};C3. 

Any symmetric configuration of the Borromean rings is the image under the pyritohedral 
group of a single embedded arc in the closed positive quadrant of the xy-plane, 
extending from a point / on the ;c-axis to a point T on the y-axis, as shown in Figure 18. 
Conversely, given any such arc IT , its images under 3*2 will form a link isotopic to 
the Borromean rings, as long as T and / are not at the same distance from the origin. 
We will assume that |/| < \T\ and will call / the intip while T is the tip. To make the 
link , the arc IT must be and must meet the axes perpendicularly at its endpoints. 

The only other points of the link in this quadrant of the -plane are 7 and f ; they will 
be important in the following discussion. 

The arcs IT of interest to us consist of a small concave piece near the intip joined to a 
large convex piece ending at the tip. That is, there is a jump point J ^ IT such that the 
arc // is strictly concave, while JT is strictly convex. As in our discussion of the clasp, 
we will parametrize // by the angle tp (less than ^ ) that its tangent vector makes to the 
right of the vertical, or by v = simp. Here v ranges from at / to some value a at J, 
which will be one of the fundamental parameters for the curves we describe. 

Along the convex arc JT we can still define v = simp, which now decreases from a 
through to — 1 . But we also use the angle if = | + ^, the angle above the horizontal 
made by the tangent vector to JT. Since our curve is , we have ip(T) = and 
ip(J) = f + arcsin a . In the curves we describe, some initial subarc JR of JT is part 
of the unit circle around 7 ; we have ^piR) < f so that along RT we can also use the 
parameter u = sirup. 

Finally, to achieve a force balance we will find it necessary that some point M along 
the circle JR has a strut to f as well as to 7. This lets us transmit some force from the 
large convex arc of one component to the smaller concave arc of another, indirectly 
through the third component. In the xy-plane, we find that M is the midpoint of the 
segment IT, and thus if we set p := sin ip{M) < cr we have 

7=(2p,0), M={p,y/\-p^), f = (0,2^1 
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Figure 18: Any configuration B of the Borromean rings with 3*2 symmetry is generated by 
a planar arc IT . We consider arcs where // is concave and JT is convex. The other points 
of B in this quadrant are the rotation images 7 and T of I and T . In our configurations, there 
are points M and R such that JR is part of the unit circle around I, and M is the midpoint 
between 7 and T . The four dotted lines are thus struts of length 1 . The height difference from / 
to 7 is cr = sin ip{J) as delineated by the horizontal dashed line, and the coordinates of M are 
given in terms of p = sin tp{M) = — cos (p(M) . 

10.2 The configuration built from circular arcs 

The configuration 62 we described in [4] is generated by an arc IT of this form. In B2 , 
we have R = T, so that the entire convex arc JT is part of the unit circle around 7. 
Furthermore, the concave arc // is also part of a unit circle, centered at f. This implies 
that M = J and a = p =: p2- The value p2 is determined by the fact that / and f 
are at unit distance, meaning 2p2 + 1 = 2\/ \ — p^. As we computed in [4], the total 
length of B2 is then Sir + 24 arcsin p2 ~ 29.0263. 

This configuration is not balanced (and thus not critical) for link-ropelength. To balance 
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the curvature forces of the circular arcs, the fans of struts to their centers would have to 
carry force proportional to arclength. But these struts would then concentrate outward 
force on the tips and inward force on the intips; there are no further struts to balance 
these forces. This is like the picture for the naive clasp — all the force is concentrated on 
the tips. As for the clasp, the tips in the critical configuration will be further apart. 

In [2], we introduce a family of thickness measures with variable stiffness. For stiffness 2 
(meaning that the curves cannot have osculating circles of diameter less than 2) we will 
see that B2 is balanced and hence critical for ropelength. Because the circular arcs have 
exactly the maximum allowed curvature, we will see that their curvature force need not 
be balanced pointwise, but only in total. Outward strut force on their midpoints (the 
tips and intips) can in a sense be spread out to balance the curvature all along the arc. 
Because p2 / 45° , however, there is an imbalance of total curvature forces between the 
convex and concave arcs. Thus our balancing measure will need an atom of force on the 
special colinear struts {/,M} and {M, T} ; this transmits force from f through M to /. 

10.3 Configurations involving clasp arcs 

To get a balanced configuration Bq of the Borromean rings, we have to replace the 
concave circular arc IJ (and part of the convex arc) by a tight clasp arc. Suppose IJ is 
part of a r-clasp for some t > a. We will now describe a configuration determined by 
certain values of our three parameters 

0<p<CJ<T<l, 

a particular curve of the class illustrated in Figure 18. 

First, the arc IJ is the piece v G [0, cr] of the r-clasp, translated out along the .x-axis 
until its tip / is at (2/5,0,0). Next, JMR is an arc of the unit circle around /, with 
v{J) = a, v(M) = p and u{R) = r. Note that to get these arcs to match up at /, we 
will need two conditions on our parameters p, a and r. Finally to define the remaining 
arc RST, consider the image IJM of IJM, rotated into the yz -plane. Then EST is 
conjugate to IJM in the sense of Lemma 9.3: it is the intersection of the unit-radius 
tube around IJM with the xy-plane, with S defined to be the point conjugate to 7. 
Figure 19 shows the arc IT and its two rotated images, that is, the part of Bq lying in 
the nonnegative orthant in space. 

Lemma 10.1 For any fixed r, suppose the parameters < p < a < t satisfy the two 
Qeometry & Topology 10 (2006) 



2106 



Cantarella, Fu, Kusner, Sullivan and Wrinkle 




Figure 19: One octant of the critical Borromean rings Bo consists of three rotated images of an 
arc IJMRT of the type shown in Figure 18. The dotted lines are struts of length 1 connecting 
the labeled points. We now describe all other struts to IT in this octant. Of course, all along the 
circular arc JMR there are struts to its center 7 . Also, between the marked struts are several 
one-parameter families of struts, joining two arcs. The first family joins the conjugate clasp 
arcs RS and 77; a second family connects 57" to the circular arc JM . The other families are 
rotated images of these, connecting JM to ST , and // to RS. The struts and {M, 7} 

are colinear. To balance IT , it is important to consider also the mirror-image struts across 
the xy-plane. This figure is an accurate drawing of Bq , except that we have exaggerated the 
separation between M and J: their actual distance is smaller than the width of the hues used in 
the picture. 



equations 

(10-1) = 2/)-Vl-cr2 + 

(10-2) = 1 - (2p - (7)2 - 



„=0 Kr(M)Vl - "2 
1-^2 



1 - a\T - af ' 

where Kr is the curvature of the clasp from Theorem 9.5. Then there is a and 
piecewise analytic arc IJMRST as described in the last paragraph. Its images under 
the symmetry group 3*2 form a configuration B{p, a, r) of the Borromean rings with 
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link-thickness LThi = 1 . 

Proof As a point on the unit circle JR around 7, the jump point / has coordinates 



Equating these, using 



from Theorem 9.5, gives (10-1) and (10-2). 

If these equations are satisfied, then the position of J is well-defined, and IJR is a 
arc, meeting the ;c-axis perpendicularly. The arc RST is the conjugate of IJM and thus 
is by Lemma 9.3. At T, the same lemma shows it meets the y-axis perpendicularly. 
At R, the u = T base of the r-clasp agrees even to second order with the unit circle. 

In this configuration, all the struts shown in Figure 19 have length 1 . If the link-thickness 
were less than 1 , there would need to be some shorter strut in this positive octant. But 
that strut would be governed by Lemma 9.3, and (rotating to assume one endpoint is 
on IT) its projection to the -plane would be normal to the arc IT; the figure makes it 
clear that no such strut exists. □ 

10.4 The balanced configuration 

Finally, we wish to find the third condition on our parameters p, a and r , which will 
ensure that B{p, a, r) satisfies the balance criterion of Corollary 6.3. 

For most of the struts, it is immediately clear what stress they need to have in a balancing 
measure p: The struts from IJ to RS, and those from MR to 7 and from RT to IM 
must be stressed exactly enough to balance the curvature force of // and of MRT. The 
conjugate clasp arcs IJ and RS exactly balance each other's curvature forces in this 
way. 

The situation along the short circular arc JM is more complicated. The struts inward 
to 7 need to balance not only the curvature force of JM itself, but also the force acting 
inward on JM from the struts from ST. Remember that the measure needed on these 
last struts is determined by the curvature of ST; this in turn determines the measure 



(a/T^, 2p-a, 0). 



As a point on the r-clasp 77, its coordinates are 
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needed on the inward struts from JM. We will write this down explicitly below. The 
final condition on our parameters then comes from a balance of forces at 7, where a 
whole family of struts converges. 

Note that this configuration Bq of the Borromean rings is the first known example of a 
ropelength-critical configuration in which this sort of transmitted force appears. Struts 
impinge on the arc JM from the direction opposite its own curvature, and transmit 
their force through that arc. Without this force transmitted through the (very) short 
arc JM, the relatively long convex piece RT would exert too much inward force on the 
relatively short concave piece //. Instead, some of this inward force, when transmitted 
through JM, becomes force outward on the concave piece 77. This transmitted force 
plays the same role in balancing Bq that the atomic force from f through M to 7 
played in balancing B2 for the stiff problem. But here our strut measure is absolutely 
continuous, with no atoms. 

To write down the final balancing condition at 7, we begin with an application of 
Lemma 9.4: the total curvature force of JMR and its mirror image across the -plane 
acts on 7 downward in the y-direction, with magnitude 

2{u{J) - u{R)) = 2{Vl -(t^-t). 

But the struts from JM carry extra transmitted force. To determine this, consider the 
curvature force of an infinitesimal arc of Sf and its mirror image across the ;cj-plane. 
Parametrizing them as usual by u. Lemma 9.4 tells us the net force, exerted in the 
negative x direction, is 2du . This horizontal force is exerted on an infinitesimal piece 
of JM and its mirror image across the -plane. If we parametrize JM by v = sin , 
then remembering that the force on this arc acts perpendicular to the arc, we see that 
if its horizontal component is du, then its vertical component is v du/Vl - v2. This 
force gets transmitted through to 7. Because of the symmetry across the yz -plane, of 
course only the vertical component matters in the end. But this symmetry also doubles 
that vertical force. (Four copies of the arc Sf act on 7: the original, and reflections 
across the xy- and -planes.) The resultant total transmitted force on 7 is upward with 
magnitude 

2 / —3l= du. 

Ju=0 

Here the upper limit of integration is u(S) = r — v(J) because / and S are conjugate 
points on the r-clasp. To make this integral explicit, we need to give the relation 
between u and v; this comes from Lemma 9.3. Along JM we have y-coordinate 2/9 — v, 
so the lemma gives 

2 , ,2 1 - (2p - vf/v ^ 

U = m(v) = 

1 - (2p - v)2 
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If one wanted, this could be solved to give v as the root of a quartic equation in p and u. 
Note that a = at v = as we expect for T and M. Plugging in a = r — a and 
V = cj (at 5 and J) reproduces (10-2). 

Summarizing, we can write the force-balancing condition at 7 as 



and so we have proved the following theorem. 

Theorem 10.2 Suppose p = po, a = ao and r = tq satisfy the three equations (10-1), 
(10-2) and (10-3). Then the configuration Bq = B{po, ao, tq) of the Borromean rings, 
constructed as in Lemma 10.1, is strongly critical for link-ropelength. 

It is easy to solve (10-1) for p, or (10-2) for p or r , or (10-3) for r, thereby eliminating 
one of our three variables. Then we are left with two nonlinear integral equations in the 
other two variables. While we have not proved formally that a solution to this system 
exists, we have solved it numerically to high precision, both in Mathematica and using 
MINPACK [26] and quadpack [30]. We obtain 



where again we follow the standard convention that the error is less than ±1 in the last 
digit shown. There is nothing delicate about this solution, since our expressions vanish 
to first order at this point. Numerically it is also clear that this solution is unique. 

Using these constants, we compute the length of our critical Borromean rings Bq 
as 29.0030. By comparison, the length of the piecewise circular Borromean rings B2 
was 29 .0263 . Thus our critical configuration Bq beats the naive circular configuration 62 
by slightly less than one-tenth of one percent. For comparison, the best lower bound 
known so far [4] for the length of the Borromean rings is 6it. 

Figure 20 shows an arclength plot of the struts in the Borromean rings. In Figure 21 we 
plot the curvature of the critical Borromean rings Bq as a function of arclength. Note 
that it is discontinuous only at J and S. Each component in Bq is built of 14 analytic 
pieces, joined in a C^'' fashion at the symmetric images of the points I, J, R and 5. 
The maximum curvature (at the intips /) is (1 — Tq)^ '^^ ~ 1.528. Therefore Bq is also 
ropelength critical for the standard ropelength functional of [4], as we will show in [2]. 
It is also critical for all the stiff ropelength functional where the lower bound A on the 
diameter of curvature is less than 2\/ 1 — Tq « 1.3. 

We note that Starostin has described [34] a configuration Bs of the Borromean rings 
with ropelength intermediate between that of our B2 and Bq ; his configuration replaces 
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I J M R ST 

Figure 20: This picture shows a portion of the strut set of our Borromean rings, plotted as pairs 
of arclength coordinates along components of the link. The horizontal axis represents arclength 
along one quadrant of the horizontal component, from I to T. On the vertical axis, we plot 
arclength along quadrants of the other two components simultaneously. (This plot accurately 
depicts the small arclength between M and /, in contrast to Figure 19 where this distance is 
exaggerated.) The horizontal segment at the bottom shows the struts from the circular arc JMR 
to 7; it joins to arcs representing the families of struts from RS to IJ and from ST to JM. 
Symmetrically, the struts to the third component are shown at the upper left: a vertical segment 
for the circle JR around /, and arcs for the struts from IJM to RST. Remembering that this 
square plot should be reflected across all of its sides to show the complete strut set, we can 
easily read off the number of stmts coming in to any point on the curve: two along // and RST, 
three along JM and one along MR . 



the arcs // and RT of B2 by clasp arcs, but does not incorporate the other features 
of Bo ■ While B5 can be balanced almost everywhere and Starostin appears to assume 
that it is a critical configuration, in fact it is not balanced at the intips since it does not 
satisfy the equivalent of (10-3). Thus by Corollary 6.3, Bs is not critical. 
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Figure 2 1 : The curvature k of the Borromean rings Bo , plotted as a function of arclength s along 
one quarter of one component of this critical configuration. The curvature has its maximum 
(about 1 .533) at the intip /, at = in this plot, and then smoothly drops off to below 1 . (This 
first part could have been plotted negatively, since this is the concave piece of Bq , but we have 
chosen to show the unsigned k of a space curve.) After a jump at 7, we have k = 1 along the 
circular arc JMR. Along the clasp arc RS, the curvature drops smoothly from 1 and then rises 
slightly again, before jumping up above 1 at 5 and then increasing to a local maximum at T . 

11 Open problems and further directions 

Our work in this paper has been motivated by a simple principle: that the ideas of rigidity 
theory for finite frameworks of bars and struts can be extended to handle mechanisms 
built from continuous curves of constraints and contacts. In the simple case of links 
critical for link-ropelength, this method has already yielded some strong results, such 
as our -regularity theorem, as well as some surprises like the tight clasp and the 
critical Borromean rings. Furthermore we expect that these methods in general, and our 
Kuhn-Tucker Theorem 5.4 in particular, will prove to be useful tools, with applications 
to a number of outstanding problems in the geometry and topology of curves and 
surfaces. 

We have mentioned our forthcoming extension of these results [2] to the classical 
ropelength problem, where the presence of curvature constraints and self-contacts of the 
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tube around individual components makes the situation considerably more challenging. 
Our theory of generalized links and obstacles should also be applicable to the study of 
packing problems for tubes and surfaces, as when thick rope is packed into a box [18] 
(a problem of some interest in molecular biology [21, 22]), or when the gray matter of 
the brain is folded and pressed against the skull. We should also mention that while 
we have only considered minimizing length in this paper, our framework should work 
equally well for other objective functionals, such as a general theory of elastic rods with 
self-contact. 

A finite-dimensional duality theorem akin to our Kuhn-Tucker theorem is one key 
step in the proof of the Unfolding Theorem of Connelly, Demaine and Rote [8]: Any 
embedded, nonconvex planar polygon admits a motion that preserves all edge lengths 
and strictly increases the distance between any two points on the polygon not already 
joined by a straight line of polygon edges. 

Our theory allows us to complete part of the proof of the (conjectured) generalization to 
smooth plane curves. Whether our methods can be made strong enough to overcome 
the formidable difficulties involved in proving a smooth unfolding theorem remains to 
be seen. 

There are several specific open questions suggested by our work above. 

Question What is the regularity of a critical curve for link-ropelength? Such curves 
are at worst and at best C''"/^ 

While we have demonstrated critical configurations of the tight clasps and Borromean 
rings, we have not attempted to prove that these configurations are minimal. 

Question Are our tight clasps and Borromean rings length-minimal in their link- 
homo topy types? 

The Euclidean-cone methods of [4] seem to hold out some hope for reducing the clasp 
problem to the case where both curves are planar, but we have not investigated this line 
of attack. 
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